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ABSTRACT:  Five  pressure  vessels  that  are  used  in  the  flight  hardware 
of  spacecraft  were  inflated  until  they  burst.  Airblast  and  fragment 
parameters  were  measured  to  evaluate  the  hazards  from  the  accidental 
rupture  of  such  tanks  during  their  handling,  checkout,  and  use.  Tank 
volumes  ranged  up  to  6  ft3,  skin  thicknesses  up  to  0.368  inches,  and 
burst  pressures  up  to  8,000  psi. 

Airblast  results  showed  that  tank  rupture  generated  shockwave 
overpressures  and  impulses  at  high  enough  levels  to  be  a  significant 
hazard  to  personnel  and  equipment.  Rupture  blast  energy  could  not  be 
equated  to  a  single  TNT  blast  yield.  The  tanks  burst  into  as  many  as 
60  fragments,  some  weighing  several  pounds.  These  fragments  attained 
maximum  velocities  of  about  1400  ft/sec. 
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Thin-walled  pressure  vessels  make  up  a  large  portion  of  the  flight 
hardware  of  missile  systems.  The  hazards  resulting  from  their 
accidental  rupture  during  tests,  checkouts,  and  use  must  be  assessed 
for  safety  considerations  and  planning  purposes.  This  report  describes 
work  done  to  determine  the  available  energy  and  potential  damage  from 
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1 .  Introduction 

1.1  Background 

Thin  walled  pressure  vessels  constitute  a  large  portion  of 
the  flight  hardware  of  missile  systems.  The  Kennedy  Space  Center  (KSC) 
is  concerned  with  the  integrity  of  such  systems  and  any  hazards  that 
result  from  their  tests,  checkouts,  and  eventual  use.  The  KSC  must 
he  able  to  determine  the  available  energy  and  potential  damage  should 
a  pressure  vessel  rupture.  Such  data  are  used  for  planning  purposes 
and  for  safety  considerations  such  as  the  protection  of  personnel  and 
facilities . 

The.  •*nergy  available  from  pressure  vessel  rupture  can  be 
determined  from  a  formula  based  on  the  isentropic  expansion  of  the 
gas  from  rupture  pressure  to  ambient  atmospheric  pressure.  A  simple 
method  to  relate  this  energy  to  damaging  shock  wave  overpressures  and 
impulses  does  not  exist j  nor  can  we  determine  the  velocity  of  tank 
fragments  and  their  energy  from  such  simple  calculations. 

•mis  report  describes  a  series  of  experiments  that  are  de¬ 
signed  to  produce  empirical  data  on  which  hazard  estimates  can  be 
based.  The  appendices  of  this  report  describe  the  beginnings  of 
theoretical  efforts  that  can  be  developed  into  methods  to  predict 
blast  parameters  and  fragment  velocities  from  rupturing  vessels. 

1 . 2  Scope  o.°  the  Experiment  and  Objectives 

Five  pressure  vessels,  or  tank3,  were  destructively  tested 
by  Inflating  them  with  nitrogen  until  they  burst.  Blast  and  fragment 
data  were  obtained.  In  addition  three  anthopomorphic  dummies  were 
placed  near  the  tanks  and  their  response  to  the  blast  and  fragments 
photographed.  Particular  information  for  each  test  tank  is  presented 
in  Table  1.1.  The  tanks  are  labeled  A  through  E  In  the  order  in  which 
they  were  fired.  They  will  be  so  identified  throughout  this  report. 

Experimental  objectives  were: 

a)  to  measure  the  blast  parameters  generated  by  the  rupture 
of  nitrogen-filled  pressure  vessels,  attempt  to  express  these  results 
in  terms  of  TNT  equivalence,  and  compare  the  results  with  the  poten¬ 
tial  energy  in  the  stored  nitrogen. 
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This  vessel  was  pressure  cycled  prior  to  rupture 

All  tanks  were  made  of  the  same  material,  titanium,  6  aluminum,  4  Vanadium  alloy. 
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b)  to  determine  the  size,  weight,  distribution,  and  initial 
velocity  of  fragments  resulting  from  tank  rupture. 

c)  to  photograph  the  response  of  anthropomorphic  dummies  to 
the  effects  of  the  bursting  pressure  vessels. 
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2 .  Experimental  Procedures 

2.1  Test  Site  and  Firing  Arena 

The  experiments  were  conducted  at  the  Naval  Ordnance  Labora¬ 
tory’s  Stump  Neck  Test  Facility  near  Indian  Head,  Maryland.  The 
experimental  setup  of  the  firing  arena,  shown  in  Figure  2.1,  was  con¬ 
structed  at  the  bottom  of  a  dry  pond.  Since  our  instrument  trailer 
was  located  only  200  feet  from  ground  zero,  the  10-foot  high  hanks 
of  the  pond  protected  it  from  fragments . 

The  setup  was  constructed  inside  a  40-foot  diameter  arena 
whose  walls  were  made  from  30  celotex  panels.  Each  panel  was  made 
from  12  sheets  of  celotex  banded  to  suitable  wood  frames  (Figure  2.2). 
A  5. 4- foot  opening  was  left  in  the  arena  wall  for  an  entrance  and 
for  the  ground  camera  station. 

A  reenforced  concrete  firing  pad  was  poured  flush  with  the 
ground  in  the  center  of  the  arena.  This  pad  was  4'  x  4'  x  1.5*.  A 
3”  pipe  led  through  the  pad  to  a  centrally  located  10"  x  10"  x  10" 
hole  in  the  top  center  of  the  pad.  High  pressure  tubing  was  fed 
through  the  3"  pipe  to  a  Junction  where  connections  to  the  test  tanks 
were  made.  The  junction  suffered  little  damage  with  this  arrangement 
and  the  same  Junction  was  used  for  the  first  four  tanks. 

Other  equipment  shown  includes  the  three  blast -gage  lines, 
the  camera  station,  the  dummies,  and  the  reflecting  screen.  Location 
and  function  of  these  equipments  will  be  discussed  later  under  the 
appropriate  heading. 

2.2  Tank  Pressurization 

2.2.1  Pressurization  System 

All  tanks  were  pressurized  with  Ng.  The  main  components  in 
the  pressurization  system,  shown  in  Figure  2.3,  are  described  below. 


The  heart  of  the  pressurization  system  was  a  one- cylinder 
Corblin  high-pres3ure  pump.  Model  A4C-600,  capable  of  pumping  approxi¬ 
mately  0.25  lbs.  of  Ng  per  minute  at  a  pressure  of  9,000  psi.  A  20- 

horsepower,  220  V.  3- phase  motor  was  used  to  drive  the  pump.  An  in¬ 
dependent  30  Kw  generator  was  used,  in  turn,  to  supply  power  for  the 
motor-pump  unit.  The  input  (supply)  pressure  to  the  pump  was  allowed 
to  vary  between  1,500  and  2,000  psi  and  was  dependent  upon  the  number 
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of  valves  opened  (tanks  used)  In  the  supply  trailer  at  any  given  time. 
Figure  2.4  shows  the  relationship  between  the  supply  pressure  and  the 
volume  of  gas  the  pump  Is  capable  of  pumping  at  a  discharge  pressure 
of  9,  OCX)  psi. 


The  N2  gas  supply  was  provided  by  a  compressed  gas  tank- 
trailer  having  a  capacity  of  35*000  ft. 3  (350  gals.).  Each  of  the  33 
tanks  in  the  trailer  emptied  into  a  common  manifold  through  individual 
valves.  Hie  discharge  from  the  manifold  was,  in  turn,  controlled  by 
a  master  valve.  Hie  average  tank  pressure  varied  between  2,200  and 
2,600  psi  and  was  dependent  upon  the  temperature  of  the  tanks.  How¬ 
ever,  at  night  when  the  tests  were  conducted  the  pressure  sometimes 
dropped  to  as  low  as  2,100  psi.  It  was  found  that  several  of  the 
individual  tank  valves  leaked,  and  as  the  pressure  in  the  tanks  asso¬ 
ciated  with  them  decreased  it  was  nearly  impossible  to  raise  the 
pressure  in  the  manifold  much  above  2,000  psi  unless  a  large  number 
of  tank  valves  were  opened  simultaneously.  Experience  showed  that  a 
manifold  pressure  of  1,500  to  2,000  psi  was  quite  satisfactory  for  the 
proper  operation  of  the  pump. 

Hie  trailer  was  connected  through  a  25 -foot  length  of 
flexible  hose  to  an  electric  solenoid  valve  at  the 
pump  site.  This  valve  was  used  to  Isolate  the  supply  trailer  from 
the  rest  of  the  system  at  tank  rupture  time  and  thus  prevented  undue 
I033  of  N2  before  the  hand  operated  valves  could  be  closed. 

The  compressed  gas  at  the  discharge  side  of  the  pump  was 

passed  through  a  water-cooled  heat  exchanger  where  the  temperature 
was  reduced  to  a  few  degrees  above  ambient.  Hie  gas  was  then  fed  to 
the  test  tank  through  approximately  180  feet  of  5/16"  ID  x  9/15"  OD 
stainless  steel  tubing,  coupled  at  20-foot  intervals.  Hie  tubing  was 
rated  at  40,000  psi  burst  pressure.  Hie  orifice  in  the  couplings  was 
7/32"  in  diameter.  Except  for  the  line  valve  located  at  the  pump 
site  these  were  the  only  constrictions  in  the  line  between  the  pump 
and  the  test  tank. 


Pressure  sensing  units  were  connected  Into  the  system  for 
each  of  the  five  tests  (Fig.  2.3).  Two  were  standard  high  pressure 
bourdon- type  dial  gages  while  the  third  was  a  bridge- type  bonded 
strain  gage  having  an  overall  accuracy  better  than  1#.  Hie  bourdon- 
type  dial  gages  had  been  calibrated  shortly  before  the  test  series 
was  begun  and  again  Immediately  afterwards^and  their  accuracy  was  of 
the  order  of  1/2#  of  full  scale.  Two  complete  sets  of  pressure  mea¬ 
suring  units  were  used  to  bracket  the  expected  burst  pressures  of  the 
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teat  tanks.  One  set  consisted  of  a  0  -  500  psl  transducer  in  conjunc¬ 
tion  with  a  0  -  1,000  psi  bourdon  dial  gage  and  was  used  with  the  two 
lower  pressure  tanks.  For  the  higher  pressure  tanks,  the  gages  were 
changed  to  ones  having  a  full-scale  capacity  of  10,000  psi.  The 
standard  dial  gage3  had  100 -pound  graduations  and  could  be  read  to 
within  50  psi  very  easily.  A  calibration  curve  for  the  10,000  psi 
line  gage  is  shown  in  Fig.  2.5.  Full  3cale  output  of  the  strain  gage 
transducers  was  about  30  mv  and  required  amplification  to  enable  them 
to  drive  a  Brush  Pen  Recorder.  The  vertical  amplifier  of  a  Tektronix 
scope  was  used  for  this  purpose.  Gage  output  was  not  recorded  for 
tanks  A  and  B  since  burst  pressure  exceeded  the  strain  gage’s  rating. 

2.2.2  Pressure  System  Operation 

The  two  low-pressure  test  tanks  (tanks  A  and  B  from  Table  1.1) 
were  pressurized  directly  from  the  supply  tank  trailer,  bypassing  the 
nump.  Pressurization  was  accomplished  by  opening  the  line  valve 
(Fig.  2.3)  3lowly  and  monitoring  the  pressure  shown  on  the  line  gages. 
Gas  was  admitted  slowly  to  insure  that  the  indicated  gage  pressure  was 
the  actual  pressure  in  the  test  tank.  The  valve  was  periodically 
closed  during  the  pressure  run  to  check  for  leaks  and  to  be  3ure  that 
the  tank  pressure  and  the  gage  pressure  were  the  same. 

At  installation,  tanks  A  and  B  were  pretested  at  200  p3l  to 
check  the  integrity  of  the  complete  system.  This  check  was  made  3  to 
4  hours  pri^r  to  the  data  run  to  allow  corrections  to  be  made. 

The  burst  pressure  of  tanks  C,  D,  and  E  was  higher  than  the 
tank- trailer  pressure.  Hence,  the  pump  wa3  used  to  raise  the  pres¬ 
sure  to  the  tank’s  bursting  point.  The  pressurization  sequence  was 
as  follows.  First,  the  test  tank  was  pressurized  to  pump-supply  - 
pressure  directly  from  the  trailer.  Meanwhile  the  pump  was  operating 
in  an  idling  mode.  Following  this  Initial  pressurization,  the  pump 
was  cut  in  and  the  tank  pressure  permitted  to  rise.  Pressure  pulsa¬ 
tions  from  the  pump  made  aocurate  reading  of  the  dial  gage  impossible. 
Therefore,  the  line  valve  was  closed  slowly  until  the  pulsations 
smoothed  out.  As  with  the  other  tanks,  the  pressure  cycle  was  inter¬ 
rupted  periodically  to  insure  that  the  tank  pressure  was  the  same  as 
that  being  read  on  the  line  gage.  This  procedure  also  allowed  the 
bonded  strain  gage  to  be  constantly  calibrated  throughout  the  pres¬ 
sure  cycle. 

The  pressure  system  was  checked  at  the  installation  of  these 
tank3  by  raising  the  pressure  to  2,000  psi.  This  was  done  3  to  4 
hours  prior  to  the  data  run. 
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FIG.  2-5  10,000  PS1  LINE  GAGE  CALIBRATION 


NOLTR  72-102 


2.2.3  History  of  the  pressurization  of  Tank3  A  and  B 

Wlien  all  systems  were  ready,  tank  A  was  slowly  pressurized 
to  460  psi  at  which  point  a  leak  developed  in  a  tank  coupling.  The 
tank  was  bled  down  and  a  new  coupling  installed.  Tank  pressure  was 
then  raised  to  500  psi  where  a  leak  developed.  The  leak  sealed  itself 
when  pressure  was  reduced  slightly.  A  third  pressure  run  was  made 
and  again  the  leak  became  evident  at  500  psi.  The  leak  was  again  self¬ 
curing  when  pressure  was  reduced.  The  tank  was  then  taken  into  the 
laboratory  and  dismantled.  We  discovered  that  the  flange  bolt3  could 
be  tightened  several  turns  after  removal  of  the  safety  wire.  Exami¬ 
nation  of  the  teflon  gasket  under  the  flange  showed  that  it  was  ex¬ 
truding  under  pressure  and  causing  the  leak. 

At  this  point,  the  tank  had  been  pressurized  to  460  psi  once 
and  to  500  psi  twice.  Yet  there  was  little  evidence  of  permanent 
deformation.  The  thick  coat  of  white  paint  that  we  had  sprayed  on 
the  tank  was  not  cracked. 

The  tank  was  reassembled  with  the  flange  bolts  snugged  down. 
The  tank  was  then  pressurized  without  incident  until  it  burst  at  625 
psi. 

Prior  to  test,  the  flange  bolt3  on  tank  B  were  snugged  down. 
This  tank  was  pressurized  to  its  burst  pressure  of  600  psi  without 
incident.  Total  time  to  pressurize  tank  B  from  0  to  burst  was  21 
minutes,  about  29  psi/minute  pressure  increase. 

2.2.4  Pre33urization  of  Tanks  C  and  D 

Tank  C  was  connected  and  raised  to  its  burst  pressure  of 
8,000  psi  without  incident.  For  this  and  the  remaining  tanks,  tank 
pressure  vs  time  was  recorded.  This  information  for  tank  C  is  plotted 
in  Figure  2.6. 

Tank  I)  was  also  inflated  to  its  burst  pressure  of  8,000  psi 
without  incident.  Tank  pressure  vs  time  during  the  inflation  cycle 
is  shown  in  Figure  2.7. 

2.2.5  pressurization  of  Tank  E 

Tank  E  was  cycled  55  times.  A  cycle  consisted  of  raising  the 
tank  pressure  from  4,000  to  5*000  psi,  holding  at  5*000  psi  for  15 
seconds,  and  then  reducing  pressure  to  4,000  psi.  The  time  required 
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FIG.  2-6  PRESSURE  VS  TIME  FOR  TANK  C 
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for  one  cycle  averaged  about  7  minutes.  At  the  end  of  the  32nd  cycle 
the  operation  was  halted  for  28  minutes  so  that  camera  film  could  be 
changed.  During  this  period,  tank  pressure  was  reduced  to  3*000  psi 
to  provide  a  greater  measure  of  safety  for  personnel  working  in  the 
vicinity  of  the  tank. 

During  the  reduction  portion  of  the  cycle,  the  Hg  was  ex¬ 
hausted  into  the  supply  trailer  to  conserve  gas. 

Following  the  55th  cycle,  tank  pressure  was  increased  to 
start  the  final  inflation  to  rupture.  At  5*200  psi,  a  leak  appeared 
and  the  pressure  was  reduced.  Hie  leak  was  self-healing  as  the  pres¬ 
sure  decreased  below  5,000  psi.  Inspection  showed  that  the  seal  on 
top  of  the  tank  was  leaking.  (Tanks  D  and  E  had  openings  at  both  top 
and  bottom.)  The  top  assembly  was  removed,  cleaned,  and  a  new  "C" 
ring  installed.  (The  MC"  ring  is  a  metal  ring  with  a  C- shaped  cross 
section.  Sealing  is  produced  when  the  ring  expands  against  two  mating 
surfaces.)  Pressure  was  again  Increased  with  a  leak  again  appearing 
at  about  5*200  psi.  Operations  were  halted  for  the  night  and  a  new 
seal  assembly  was  fabricated  the  following  day.  This  assembly  U3ed 
the  more  conventional  double  neoprene  "O”  ring  seals. 

The  test  was  resumed  the  following  evening  and  no  leaks  ap¬ 
peared.  However,  at  5*250  psi,  a  signal  light  indicated  that  the 
breakwire  around  the  tank  had  parted.  Hie  pressure  was  reduced  to 
3,000  psi  and  the  breakwire  tested.  The  breakwire  was  still  intact 
and  we  assumed  that  a  voltage  transient  had  activated  the  signal 
light.  Inflation  was  resumed  and  the  burst  pressure  of  8,130  psi  was 
reached  without  further  Incident.  The  final  inflation  pressure  vs 
time  for  tank  E  is  shown  in  Figure  2.8. 

2.3  Instrumentation 

2.3.1  Blast  Instrumentation 

The  blast  instrumentation  system  is  built  around  piezoelectric 
gages  whose  outputs  are  recorded  by  cathode-ray  oscilloscopes  and  ro¬ 
tating  drum  cameras.  The  oscilloscopes  and  auxiliary  equipment  used 
for  time  bases,  synchronization,  and  power  are  housed  In  a  special 
Instrumentation  trailer.  This  and  similar  systems  have  been  used  at 
N0L  over  the  past  several  years  and  are  described  in  reference  (1)7 

Both  side- on  and  face- on  pressures  were  measured  using  three 
gage  types.  The  side-on  gages  were  either  Atlantic  Research  Corpora¬ 
tions  model  LC-33  gages  or  Susquehanna  Instruments  model  ST-7  gages. 

The  face-on  gages  (Atlantic  Research  Corporation  Model  LC-71)  were 
mounted  in  18-inch  round  aluminum  baffles  oriented  normal  to  the 
direction  of  the  shock  flow. 


y~  -  ■  — . . . -  ■ 

References  are  found  on  page  103 
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FIG.  2-8  PRESSURE  VS  TIME  FOR  TANK  E 
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. .  .  Signals  frcm  the  LC-33  and  LC-71  gages  were  transmitted  via 

anti-noise  coaxial  cable  to  charge  amplifiers  mounted  in  the  instrument 
trailer.  Signals  from  the  ST-7  gages  were  fed  to  charge  amplifiers 
located  at  the  gages.  The  voltage  outputs  from  all  charge  amplifiers 
on  cathode-ray  oscilloscopes.  Rotating  drum  cameras 
with  a  suitable  time  base  produced  permanent  pressure- time  records. 

a5e™-£recjuenc5r  resP°nse  the  system  was  from  below  0.2  hz  to  about 
oO  Khz . 

In  operation,  the  equipment  was  setup  and  ready  with  cameras 
open  during  the  inflation  cycle.  The  oscilloscopes  were  unblanked  at 
tank  rupture  by  a  signal  from  a  trigger  unit  initiated  by  a  breakwlre 
around  the  tank. 


2.3.2 


lent  Velocity  System 


^  Attempts  to  measure  fragment  velocities  involved  two  techniques. 
One  method  used  breakwires  and  electronic  counters.  The  first  break- 
wire  was  wrapped  around  the  tank  under  test  and  signaled  tank  rupture. 

The  second  breakwlre  was  supported  by  an  8~gage  wire  frame  formed  to 
the  tank  shape  but  everywhere  nominally  1  ft  larger  in  radius.  At 
tank  rupture,  the  first  breakwlre  interrupted  a  current  to  signal  an 
electronic  counter  to  start.  The  first  fragment  to  travel  the  1  ft 
separation  broke  the  second  wire.  Thi3  again  interrupted  a  current 
to  signal  the  counter  to  stop.  To  avoid  errors  in  the  counters  due  to 
trigger  level  settings,  trigger  signal  shapers  were  used  between  both 
the  start  wire  and  stop  wire  signals  and  the  appropriate  counter  input. 


The  second  technique  used  to  measure  fragment  velocities  in¬ 
volved  a  stroboscopic  photographic  system.  Briefly,  this  system  used 
two  camera  stations.  Each  station  consisted  of  a  single  4”  x  5" 
camera  and  an  electronic  flash  tube.  The  tube  was  flashed  four  times 
at  measured  intervals.  The  first  flash  occurred  at  tank  rupture, 
with  the  three  succeeding  pulses  occurring  at  later,  appropriately 
spaced  intervals.  The  strobe  unit  is  diagrammed  in  Figure  2.9. 


One  of  the  cameras  was  placed  at  ground  level  and  photographed 
against  the  background  formed  by  a  reflecting  screen.  Limiting  the 
size  of  the  light  source  in  conjunction  with  the  above  screen  produced 
a  true  shadow-graph  system  similar  to  that  described  in  reference  (2). 
Thus  In  most  cases,  the  shock  generated  by  the  bursting  tank  was  photo¬ 
graphed  by  one  or  more  of  the  flashes .  With  this  system,  however, 
only  the  shadow  of  the  fragment  against  the  screen  was  photographed. 


A  second  camera  was  located  on  a  tower  above  and  off  to  the 
the  talk  .  This  camera  photographed  tank  fragments 

directly  against  a  low  reflecting  background. 


equipment  tests,  the  flash  tubes  fired  at  the  preset  times 
within  ±  0.3ft  To  be  sure  that  they  operated  properly  at  tank  rupture. 
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FIG.  2-9  STROBE  UNIT 
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we  monitored  the  light  output  by  placing  a  photocell  near  each  flash 
tube.  Hie  photocell  outputs  were  recorded  on  a  dual  beam  oscilloscope 
with  a  suitable  time  base. 

Hie  system  was  operated  as  follows.  When  the  internal  tank 
pressure  reached  the  specified  working  level  during  the  inflation 
cycle,  the  camera  shutters  were  opened  remotely.  At  tank  rupture, 
the  first  breakwire  signal  triggered  the  flasher  system  3hown  in  Fig¬ 
ure  2.9.  At  the  same  time  the  camera  shutter  closing  solenoids  were 
energized,  and  the  shutters  closed  18  milliseconds  after  tank  rupture. 
Hie  time  from  tank  rupture  to  the  last  strobe  flash  was  always  less 
than  18  milliseconds. 

Since  the  camera  shutters  remained  open  for  several  minutes, 
the  experiments  were  done  on  moonless  nights. 

Both  systems  did  not  function  on  every  tank  rupture.  However, 
one  or  the  other  system  always  operated  so  that  fragment  velocities 
were  measured  on  every  test. 


2.3.3 


the  Anthropomorphic  Dummies 


Hie  reactions  of  the  anthropomorphic  dummies  to  tank  rupture 
were  photographed  using  16- mm  Mitchell  cameras  operating  at  24  frames/ 
second.  Since  the  experiments  were  done  at  night,  light  was  supplied 
by  photoflood  lamps .  These  lamps  were  turned  on  at  tank  rupture  to 
avoid  fogging  the  strobe  camera  film.  To  shorten  turn-on  time,  the 
filaments  of  the  photoflood  lamps  were  preheated  to  a  dull  red  glow. 
Hie  photofloods  reached  90#  of  full  brightness  in  about  50  milli¬ 
seconds  after  tank  rupture  --  20  milliseconds  relay  closure-time  and 
30  milliseconds  for  the  lamp  filaments  to  reach  Incandescence.  Thus, 
the  dummies  were  photographed  within  two  frames  after  tank  rupture . 


follows : 


2.3.4  Sequence  of  Operations 

A  typical  sequence  of  events  for  instrumenting  tank  rupture 


a)  When  the  pressure  in  the  tank  under  test  reached  its 
operating  pressure,  the  following  operations  were  done  manually: 

1)  Open  strobe  camera  shutters 

2)  Check  ready  position  of  all  electronics 

3)  Open  drum  camera  shutters  on  blast  record¬ 
ing  equipment 

4)  Start  Mitchell  cameras 

b)  At  tank  rupture,  the  following  operations  were  initiated 
automatically  by  the  first  breakwire  signal: 

1)  Unblank  oscilloscopes  in  pressure  instrumentation 
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2)  Initiate  fragment  velocity  screen  to  start 

!  fragment  velocity  counter. 

Trigger  first  strobe  flash. 

Turn  on  photofloods. 

Energize  strobe  camera  solenoid  for  shutter  closure 

c)  After  tank  rupture  but  less  than  0+1  millisecond, 
fragment  velocity  screen  pulse  stops  fragment  velocity  counter . 

d)  At  0  +  1.5  milliseconds,  second  strobe  pulse  triggers. 

e)  At  0  +  3.0  milliseconds,  third  strobe  pulse  triggers. 

f )  At  0  +  5.0  milliseconds,  fourth  and  last  strobe  pulse 

triggers . 

g)  At  0  +  18  milliseconds,  strobe  camera  shutters  close. 

h)  At  0  +  50  milliseconds,  photoflood  lamps  reach  90#  of 

maximum  brightness . 

i)  Equipment  turned  off  manually  at  about  0  +  1  minute. 
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3 .  Blast  Measurements 

3.1  Background  and  Theory 

The  purpose  of  the  blast  measurements  was  to  determine  an 
empirical  explosion  energy  output;  hopefully,  a  single  TNT  energy 
equivalent  to  aid  in  estimating  the  damage  potential  from  3uch  acci¬ 
dental  tank  ruptures. 

The  expected  side-on  airblast  overpressures  were  calculated 
to  aid  in  experimental  planning,  lhe  calculations  used  the  one¬ 
dimensional  hydrocode  from  reference  (3).  Results  of  these  calcula¬ 
tions  are  given  in  Appendix  A.  The  hydrocode  used,  called  WUNDY, 
gives  the  results  for  a  free  air  situation.  Since  the  tanks  were  rup¬ 
tured  near  the  ground  all  measurements  were  made  in  the  maeh  region. 
Therefore,  a  reflection  coefficient  of  1.5  was  used  when  comparing 
theory  and  measurements.  This  will  be  discussed  further  in  the  pre¬ 
sentation  of  the  data. 

Nominal  8.2-lb  TNT  spheres  were  fired  to  validate  our  TNT 
predictions  and  to  generate  data  to  directly  compare  with  the  data 
generated  by  the  rupture  of  the  high  pressure  helium  tanks. 

All  the  tanks  were  placed  with  their  center  of  mass  16" 
above  the  ground.  The  TNT  spheres  were  fired  with  their  center  of 
mass  12”  above  the  ground.  The  TNT  spheres  were  fired  at  thi3  level 
to  assure  that  the  blast  measurements  were  made  in  the  mach  region. 

3.2  Blast  Page  Layout 

The  blast  gages  were  located  along  three  equally  spaced 
radial  lines  extending  from  Ground  Zero  (GZ).  This  set  up  is  shown 
in  Figure  3.1.  Gage  distances  are  those  used  for  tanks  D  and  E.  A 
similar  set  up  was  used  for  tsmks  A,  B,  and  C  except  that  the  dis¬ 
tances  from  GZ  to  the  gages  were  reduced.  Nomimal  gage  distances 
for  tanks  A,  B,  and  C  were  3.5*  5.6,  and  12.0  ft  from  GZ.  All  gages 
were  6  Inches  above  the  ground. 

3.3  Pressure- Time  Histories 


Tracings  of  the  pressure-time  histories  are  shown  in  Figures 
3.2  through  3.6  for  tanks  A  through  E,  respectively.  The  traces  are 
similar  to  those  from  the  TNT  charge  shown  in  Figure  3.7.  There  sure 
incidences  of  double  peaks  and  other  perturbations  on  tank-rupture 
records,  particularly  from  the  gages  nearest  GZ.  Records  from  the 
positions  furthest  from  the  charge  are  usually  free  from  such  pertur¬ 
bations  indicating  that  the  shock  wave  has  smoothed  by  the  cirne  it  has 
propagated  this  distance. 

Shockwave  arrival  times  are  different  from  the  different  gage 
lines.  This  shows  a  lack  of  symmetry  for  the  shock  front  that  is  fur¬ 
ther  borne  out  in  the  pressure  differences  shown  in  section  3.4. 
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Gages  at  the  18  ft  positions  (tanks  D  and  E)  were  two  feet 
nearer  the  charge  than  the  celotex  arena  wall.  Therefore,  there  is 
a  strong  reflection  from  the  wall.  This  reflection  is  the  second 
shock  shown  on  the  records  from  Tanks  D  and  E  for  the  18  ft  position. 
(Figures  3.5  and  3.6). 

3. 4  Blast  Results 

3.4.1  Pressure  Results  from  Tanks  A  and  B 

Blast  data  generated  by  the  rupture  of  Tanks  A  and  B  are 
given  in  Tables  3.1  and  3*2.  Side-on  pressure  data  are  presented  in 
Figures  3*8  and  3.9*  The  predicted  TNT  curve  was  obtained  from  the 
calculated  data  in  reference  (4)  by  adjusting  for  ground  reflection. 
Specifically,  a  ground  reflection  coefficient  of  1.5  was  assumed. 

The  reference  (4)  data  were  scaled  to  0.138  pounds  TNT  by  reducing 

distances  by  (1.5)1//3  (0,138)1//3  or  0.592.  These  data  are  shown  in 
Figure  3-8.  The  measured  TNT  data  was  scaled  from  the  8.2  lb  TNT  data 
shown  in  Table  3.6.  The  figure  0.138  pounds  of  TNT  comes  from  a 
calculation  based  on  the  isentropic  expansion  of  the  gas  in  the  tank 
from  burst  pressure  to  ambient  (14.7  psi).  The  detonation  energy 
of  TNT  was  assumed  to  be  1018  calories/gram.  A  similar  isentropic 
calculation  for  the  Tank  B  rupture  pressure  and  volume  yields  a  TNT 
equivalent  of  0.169  pounds  of  TNT.  The  data  for  Tank  B  are  shown 
in  Figure  3-9.  The  predicted  TNT  curve  in  Figure  3-9  was  generated 
by  multiplying  the  distances  in  the  reference  (4)  data  by 

(1.5)1/3  (0.169)1/3  or  0.633. 

The  predictions  for  the  tank  rupture  pressures  were  obtained 
from  the  calculated  values  shown  in  Appendix  A  by  appropriate  allowance 
for  ground  reflection.  Again,  a  reflection  coefficient  of  1.5  was 
assumed.  Thus  distances  for  a  given  pressure  were  multiplied  by 

(1.5) 1/3  or  1.145. 

Two  things  are  immediately  apparent  from  Figures  3.8  and  3-9. 

The  pressure  decay  slopes  are  different  for  the  curve  based  on  tank- 
rupture  pressures  and  the  TNT  based  curve.  Secondly,  pressures  along 
different  radii  from  GZ  are  not  equal  showing  a  lack  of  symmetry. 
Therefore,  a  single  value  for  the  TNT  equivalency  of  the  tank  rupture 
energy  does  not  describe  the  phenomena. 

The  predicted  curve  was  based  on  a  burst  pressure  of  460  psig; 
inserting  the  correct  burst  pressure  in  the  calculations  will  increase 
the  predicted  side-on  pressures  from  tank  rupture  by  about  15%  in  the 
pressure  range  of  interest.  Further,  no  allowance  was  made  in  tno  cal¬ 
culations  for  the  reduction  of  energy  available  to  the  shock  resulting 
from  the  energy  expended  in  tank  rupture  and  fragment  acceleration. 
Including  this  factor  for  conventional  explosions  results  in  a  reduc¬ 
tion  of  side-on  pressure  in  the  high  pressure  region.  The  inclusion 
of  this  "Mass  Effect"  in  the  tank  calculations  may  give  similar  results. 
The  quantitative  effect  is  not  known  in  the  absence  of  such  calcula¬ 
tions  (Ref.  (5)). 
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FIG.  3-8  SIDE-ON  OVERPRESSURES  FROM  TANK  A;  1 .34  FT3  BURST  AT  625  PSI 
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Pace-on  pressures  are  shown  in  Figures  3.10  and  3.11  for 
Tanks  A  and  B.  The  TNT  curve  was  generated  from  reference  (4)  data 
using  the  same  scale  factoisas  for  the  side-on  pressures.  Pressures 
were  calculated  using  the  equation  for  normally  reflected  pressures 
given  in  reference  (6). 

3.4.2  Pressure  Results  from  Tank  C 

Data  from  this  tank  are  presented  in  Table  3.3  and  Figures 
3.12  and  3.13.  The  predicted  TNT  curve  was  based  on  a  yield  of  0.324 
lbs  and  a  reflection  coefficient  of  1.5.  Thus  the  TNT  curve  was  ob¬ 
tained  by  decreasing  the  distances  for  the  calculated  one-pound  data 

given  in  reference  (4)  by  (l.S)1^  x  (Q.324)1^  or  0.79.  No  TNT 
charges  of  equivalent  size  were  fired  since  the  8.2-lb  TNT  data  found 
in  Table  3.6  can  be  scaled  down  for  comparison.  The  predicted  curve 

for  tank  rupture  pressures  was  obtained  from  the  6~ft^  vessel  pre¬ 
dictions  in  Appendix  A  as  follows.  Since  the  burst  pressure  of  the 
tanks  were  equal,  distances  for  equal  pressures  were  reduced  by  the 

4  6  . 

ratio  of  the  tank  radii  — —  or  0.34.  The  tank  rupture  pressure  pre- 

13.5 

dictions  did  not  allow  for  the  energy  losses  due  to  tank  rupture  and 
acceleration.  An  accounting  for  this  energy  in  the  calculations  will 
lower  side-on  pressures  nearest  the  tank;  the  amount  is  unknown  until 
such  calculations  are  made. 

The  measured  bank  rupture  data  again  show  a  lack  of  symmetry 
about  GZ  and  a  different  slope  from  that  of  the  predicted  TNT  curve* 

Thus  a  single  TNT  energy  equivalent  does  not  adequately  describe  the  data. 

Similar  comments  are  appropriate  to  the  face-on  pressures 
shown  in  Figure  3.13.  The  differences  between  measured  data  and  pre¬ 
dicted  TNT  data  are  accentuated.  This  comes  about  because  the  ratio 
between  3ide-on  and  reflected  pressures  is  pressure  dependent.  It 
varies  from  a  low  of  2  as  side- on  pressures  approach  zero  to  8  as  they 
become  infinitely  large. 

3.4.3  Pressure  Results  from  Tanks  D  and  E 

Pressure  data  from  the  rupture  of  these  two  tanks  are  pre¬ 
sented  in  Tables  3.4  and  3.5  and  Figures  3.14  and  3.15. 

The  predicted  TNT  curve  in  Figure  3.14  is  based  on  the  1-lb 
free  air  data  in  reference  (4).  A  ground  reflection  coefficient  of 
1,5  was  assumed.  Thus  the  predictions  for  10.1  lbs  of  TNT  were 
obtained  by  multiplying  the  distances  for  1  lb  by  the  scale  factor 
(1.5) Vs  (10.1)1/*  or  2.48.  (The  tank  rupture  energy  based  on  lsen- 
troplc  expansion  is  equivalent  to  10.1  lbs  of  TNT.)  Measured  8.2-lb 
TNT  data  were  scaled  to  10.1  lb3  by  increasing  distances ( 10. 1\1/3 
or  1.07.  ylfiP 

The  predicted  tank  rupture- pressures  were  obtained  by  mul¬ 
tiplying  the  Wundy  code  calculations  by  a  suitable  factor  to  account 
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FIG.  3-10  FACE-ON  PRESSURES  FROM  TANK  A;  1 .34  Ft3  BURST  AT  625  PSI 
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FIG.  3-11  FACE-ON  PRESSURES  FROM  TANK  B;  1 .68  FT3  BURST  AT  600  PSI 
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FIG.  3-12  SIDE-ON  OVERPRESSURES  FROM  TANK  C;  0.235  FT3  BURST  AT  8,000  PSI 
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FIG.  3-14  SIDE-ON  OVERPRESSURES  FROM  TANKS  D  AND  E;  6  FT3 
BURST  AT  8,000  PSI 
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FIG .  3-15  FACE-ON  PRESSURES  FROM  TANK  D;  6  FT3  TANK  BURST  AT  8,000  PSI 
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for  pressure  enhancement  due  to  the  ground  proximity.  As  with  TNT,  a 
reflection  coefficient  of  1.5  was  assumed.  Therefore  the  scale  factor 

is  (1.5)  or  1.145.  Again,  the  effects  of  the  pressure  vessel's 
mass  was  not  considered.  Inclusion  of  the  mass  effect  in  the  side-on 
pressure  calculations  for  tank  rupture  would  give  results  similar  to 
those  predicted  for  tanks  A,  B,  and  C. 

Figure  3.14  shows  the  measured  tank  rupture  pressures  from 
tank  D  to  he  in  reasonable  agreement  with  those  from  tank  E.  Again, 
the  shock  strength  is  not  symmetrical  about  GZ.  There  is  also  a  dif¬ 
ference  between  the  slopes  of  the  predicted  TNT  curve  and  that  ob¬ 
tained  from  the  tank  rupture  data."  Thus,  no  single  TNT  equivalence 
value  can  be  used  to  describe  the  tank  rupture  data. 

Similar  comments  apply  to  the  face-on  pressures  shown  in 
Figure  3.15«  Hie  differences  between  predicted  face-on  and  measured 
face-on  pressures  are  greater  than  those  for  the  side-on  case.  Again, 
this  comes  about  because  the  ratio  of  face-on  to  side-on  pressures  is 
pressure  dependent. 

3.5  Positive  Shockwave  Durations 

3.5.1  Durations  for  Tanks  A  and  B 

Figures  3.16  and  3.17  present  the  positive  shock  wave  dura¬ 
tions  from  Tanks  A  and  B.  Predicted  TNT  values  were  scaled  from 
reference  (6)  data  to  the  isentroplc  expansion  energy  eauivalents  of 
0.138  lbs  and  O.169  lbs  TNT  for  Tanks  A  and  B  respectively.  Both 
distances  and  durations  were  multiplied  by  O.592  for  Tank  A  and  O.633 
for  Tank  B.  The  origin  of  these  scale  factors  is  given  in  Section 
3.4.1.  The  measured  8.2-lb  TNT  data  were  scaled  for  Tank  A  by  mul¬ 
tiplying  measured  distances  and  durations  by  the  factor^O.138^1/* 

or  O.256  and  to  Tank  B  by  the  factor  /0.l69v  V *  or  0.275?*2 

v— ) 

Figures  3.16  and  3.17  show  that  tank  rupture  durations  are 
longer  than  those  from  TNT.  Because  the  tank  ruoture  durations  are 
much  longer  close-in,  they  generate  a  curve  with"  a  shallower  slope  than 
does  TNT. 

3.5.2  Durations  for  Tank  C 

Figure  3.18  presents  the  positive  shock  wave  durations  gene¬ 
rated  by  the  rupture  of  ^ank  C.  Predicted  TNT  values  are  those  from 
reference  (6)  scaled  to  0.324  lbs  TNT  burst  near  a  reflecting  surface 
with  1.5  reflection  coefficient.  Thus  both  distances  and  durations 
were  multiplied  by  (l.5)V  (0.324)1/3  or  O.79.  Measured  8.2-lb  TNT 

data  were  scaled  to  0.324  lbs  by  multiplying  both  durations  and  dis¬ 
tances  by,0.325\V3  or  0.34. 
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FIG.  3-17  SIDE-ON  POSITIVE  DURATION  FROM  TANK  B;  1 .68  FT3 
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Figure  3.18  shows  longer  positive  shock  wave  durations  than 
those  predicted  for  0.324  lbs  of  TNT  at  the  closer  positions.  Thus 
the  tank  rupture  data  generates  a  duration  curve  that  has  a  flatter 
slope  than  does  the  curve  generated  by  a  conventional  explosion. 

3.5.3  Durations  for  Tanks  D  and  E 

Figures  3.19  and  3.20  present  the  side-on  positive  shock 
wave  durations  from  the  larger  helium  vessel  rupture,  ihe  predicted 
TNT  curve  is  generated  from  reference  (6).  These  data  were  scaled  to 
10.1 Its  of  TNT  fired  on  the  ground  by  assuming  a  ground  reflection 
coefficient  of  1.5.  The  experimental  TNT  values  are  those  measured 
from  the  8.2-lb  TNT  charges  scaled  to  10.1  lbs. 

At  the  18  ft  gage  position,  the  wall-reflected  shock  arrived 
back  at  the  gages  before  the  incident  wave  had  decayed  to  zero  (Figure 
3.5).  Therefore,  positive  durations  for  this  position  were  obtained 
by  extrapolating  the  incident  pressure  history  to  zero  pressure. 

The  pattern  of  longer  positive  shock  wave  durations  at  the 
close-in  gage  positions  was  present  in  these  two  events  as  it  was  in 
all  tank  rupture  cases. 

3.6  Positive  Shockwave  Impulse 

3.6.1  Impulse  from  Tanks  A  and  B 

Figures  3.21  and  3.22  present  the  positive  shock  wave  impulse 
generated  by  the  rupture  of  tanks  A  and  B.  The  predicted  TNT  data 
are  scaled  from  the  reference  (6)  data  using  an  assumed  reflection 
coefficient  of  1.5.  Thus  the  one- lb  TNT  data. from  reference  (6),  both 
impulse  and  distance,  are  multiplied  by  0.592  or  0.633,  the  scale 
factors  developed  in  Section  3.4.1  for  Tanks  A  and  B  respectively. 

The  measured  8.2-lb  TNT  data  were  sealed  to  these  tanks  also.  This  was 
done  by  multiplying  the  measured  distances  and  impulses  by  0.256  or 
0.275,  the  scale  factors  given  In  Section  3.5.1  for  scaling  the  8. 2 -lb 
data  to  Tanks  A  and  B  respectively. 

The  tank  rupture  data  are  higher  than  the  predicted  TNT  data 
and  except  for  the  farther  out  position,  higher  than  the  measured  TNT 
data. 

Face-on  positive  impulses  have  not  been  plotted  but  are  given 
in  Tables  3.1,  3.2,  and  3.6  for  tanks  A,  B,  and  TNT  respectively. 

3.6.2  Impulse  from  Tank  C 

Figure  3.23  presents  the  positive  shock  wave  impulse  data 
from  the  rupture  of  tank  C.  The  predicted  TNT  curve  was  scaled  to 
0.324  lb3  from  the  data  in  reference  (6)  using  an  assumed  reflection 
coefficient  of  1.5.  Thus  the  referenced  distances  and  impulses  were 

reduced  by  a  factor  of  il.5)1^  (0.324)1/3  or  0.79.  The  8. 2 -lb  TNT 
data  were  scaled  to  0.324  lbs  by  reducing  both  distances  and  impulses 

or  0.34. 
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FIG.  3-19  SIDE-ON  POSITIVE  DURATION  FROM  TANK  D;  6  FT3 
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FIG.  3-20  SIDE-ON  POSITIVE  DURATION  FROM  TANK  E;  6  FT3 
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FIG.  3-21  SIDE-ON  POSITIVE  IMPULSE  FROM  TANK  A;  1 .34  FT3 
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FIG.  3-22  SIDE-ON  POSITIVE  IMPULSE  FROM  TANK  B;  1 .68  FT3 
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FIG.  3-23  SIDE-ON  POSITIVE  IMPULSE  FROM  TANK  C;  0.235  FT3 
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Side-on  positive  impulse  from  tank  rupture  seems  to  scatter 
about  the  predicted  TNT  data.  Pace-on  positive  impulse  was  not  plotted 
but  is  given  in  Table  3.3. 

3.6.3  Impulse  from  Tanks  D  and  E 

Figures  3.24  and  3.25  present  the  positive  shock  wave  impulse 
data  from  the  rupture  of  tanks  D  and  E.  The  predicted  impulse  for  TNT 
is  generated  from  the  1-lb  reference  (6)  data  by  assuming  a  reflection 
coefficient  of  1.5.  Thus  the  referenced  data,  both  distance  and  impulse, 

were  increased  by  (1.5)  V  3  (10.1  )V*  or  2.48.  The  measured  6.2-lb 

data  scaled  to  10.1  lbs  is  seen  to  compare  quite  well  with  the  predicted 

curve 

The  tank  rupture  data  generally  scatter  about  the  predicted 
TNT  curve.  However,  the  impulse  measured  along  line  1  for  tank  D 
(Figure  3.24)  are  quite  low.  The  side- on  pressures  along  this  line 
for  this  tank.  Figure  3.14,  are  also  quite  low. 
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FIG.  3-24  SIDE-ON  POSITIVE  IMPULSE  FROM  TANK  D;  6  FT3 
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FIG.  3-25  SIDE-ON  POSITIVE  iMPULSE  FROM  TANK  E;  6  FTJ 
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4.  Fragmentation  Results 

4.1  Sat  up  for  Fragment  Samples 

An  arena  40  ft  in  diameter  was  used  to  catch  a  sampling  of 
the  fragments.  Arena  walls  were  constructed  of  panels  made  from 
layers  of  4«  x  8*  x  celotex  wallboard. 

The  thickness  of  the  panels  (number  of  pieces  in  each  wall) 
was  calculated  to  stop  most  fragments.  Prom  reference  (7),  depth 
of  fragment  penetration  can  be  calculated  from  the  following  empirical 
relationship: 

.23 


where  dp  =  penetration  depth  in  inches  normal  to  panel  surface 
v  *  fragment  velocity  in  ft/sec 
M  «=  weight  of  fragment  in  grains 

2 

A  =  presented  area  of  fragment  in  in 

0  =  angle  between  fragment  trajectory  and  a  normal  to  the 

panel  surface 

Not  having  factual  fragmentation  data  available  from  rupturing  pressure 
vessels,  for  our  case  we  considered  a  fragment  that  is  1"  x  1"  x  0.3" 
striking  the  panel  in  such  a  manner  as  to  present  an  area  of  about 
0.5  in  .  The  fragment  weighs  330  grains  and  strikes  the  panel  at  a 
velocity  of  1000  ft/sec.  Since  the  greatest^panel  thickness  is 
required  when  the  trajectory  is  normal  to  the  panel  surface,  sec  9 
was  taken  as  1.0.  For  the  above  conditions,  the  penetration  depth 
is  about  5.0".  Based  on  this,  panels  6"  thick  should  be  adequate  to 
catch  most  fragments. 

The  total  number  of  panels  used  for  the  40  ft  arena  was 

=  10^  =  31.4  panels.  In  practice,  cnly  30  panels  were  used  with  the 

space  equivalent  to  1.4  panels  (5.6*)  left  open  for  an  entrance  and 
for  a  camera  port.  The  thirty  panels  that  made  up  the  arena  wall  were 
numbered  1  through  30  starting  at  the  arena  opening  and  proceeding 
in  a  clockwise  direction.  The  total  number  of  celotex  sheets  used  was 
12/panel  for  30  panels  or  360  sheets.  Details  of  arena  construction 
are  shown  in  Figure  2,2. 

To  complete  the  tank  experiments  in  a  rather  restricted  time 
frame,  construction  materials  were  ordered  on  the  above  preliminary 
information.  Later  and  more  accurate  information  proved  that  fragment 
material  density  was  greater  than  that  used  in  the  calculations. 

Thus  while  the  arena  walls  proved  capable  of  stopping  fragments 
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from  the  thin-walled  cylindrical  tanks,  many  of  the  fragments  from 
the  higher  pressure  tanks  penetrated  the  arena  walls.  S<jme  of  those 
that  penetrated  the  walls  were  recovered. 

4.2  Fragment  Recovery 

4.2.1  Fragments  from  Tank  A 

A  total  of  36  fragments  of  tank  A  were  recovered.  These 
are  shown  photographed  over  a  calibration  grid  in  Figure  4.1.  Grid 
line  spacing  in  this  and  all  fragment  photographs  is  2  inches. 

Fragment  sizes  range  from  about  10"  x  12"  to  1"  x  2".  Weights  for  each 
fragment  are  given  in  Table  4.1.  Those  fragments  that  were  found 
imbedded  in  the  arena  walls  are  identified  by  the  panel  number  in  which 
they  were  found  in  Table  4.1. 

Table  4.2  lists  the  number  of  hits  cn  panels  for  each 
shot.  A  check  shows  that  the  number  of  panel  hits  for  shot  A  does  not 
correspond  to  the  10  fragments  found  in  the  celotex  panels  as  listed 
in  Table  4.1,  The  reason  for  this  is  that  Table  4.2  lists  the  total 
panel  hits,  those  fragments  that  remain  in  the  panel  and  those  that 
rebound;  Table  4.1  lists  only  those  fragments  that  remain  in  the 
celotex. 

The  remaining  fragments  listed  in  Tab3e  4.1,  except 
for  23,  were  found  scattered  about  the  arena.  Fragment  #23  was  found 
150  ft  from  GZ  in  a  radial  line  with  panel  6. 

Prior  to  firing,  the  tank  weighed  3850  grams.  Total 
weight  of  the  fragments  found  was  2837  gf^ms,  about  73 %  of  the  total 
weight.  This  recovery  percentage  may  distort  the  fragment  picture 
since  the  recovered  items  include  the  heavy  tank  parts  listed  in 
Table  4.1  and  Figure  4. 1,  fragments  1,  34,  35*  and  36.  The  amount  of 
tank  skin  recovered  weighed  about  1349  grams.  The  total  weight  of 
the  tank’s  skin  based  on  a  density  of  4.41  and  a  thickness  of  0.02 
Inches  was  2390  grams.  Thus  about  56$  of  the  tank's  skin  was  recovered. 

4.2.2  Fragments  from  Tank  B 


A  total  of  36  fragments  from  Tank  B  were  recovered.  These 
are  shown  in  Figure  4.?.  Sizes  range  from  about  16"  x  17"  down  to 
about  2"  x  l",  (Fragments  2  and  33.)  Table  4.3  gives  the  weight  of 
each  fragment.  Fragment  numbers  correspond  to  those  in  Figure  4.2. 

Table  4.2  indicates  that  a  total  of  15  fragments  from 
Tank  B  hit  the  arena  walls;  of  these,  10  were  identified  and  are  listed 
in  Table  4.3. 


The  remainder  of  the  fragments,  except  for  #4,  were  found 
scattered  about  the  arena  floor.  Fragment  #4  was  found  about  100  ft 
from  GZ  and  in  line  with  panel  26. 
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FIG.  4-2  FRAGMENTS  FROM  TANK  B  (2"  LINE  SPACING) 
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Prior  to  firing  the  tank  weighed  4620  grams;  total 
weight  of  the  fragments  found  was  3300  grains— about  71#  of  the  tank 
weight,  out  of  a  calculated  total  skin  weight  of  2890  grams,  we 
recovered  177 ^  grams  or  6l#, 

4.2.3  Fragments  from  Tank  C 

The  8  fragments  recovered  from  Tank  C  are  shown  in 
Figure  4.3.  Fragnent  sizes  range  from  about  4M  x  6M  down  to  about 
1”  x  2".  Weights  for  each  fragment  are  shown  in  Table  4.4. 

Three  panel  hits  are  shown  for  Tank  C  in  Table  4.2. 
These  three  fragments  were  located  and  are  listed  in  Table  4.4. 

The  total  weight  of  fragments  recovered  for  this  tank 
was  617  grams.  The  tank  weight  at  firing  was  2860  grams.  Thus  only 
21.5#  of  the  total  weight  was  recovered.  The  reason  for  the  low 
recovery  percentage  may  have  been  due  to  the  ricochet  of  many  of  the 
fragments  off  the  concrete  pad  on  a  trajectory  that  carried  them  over 
the  arena  wall. 

4.2.4  Fragments  from  Tank  D 

Figure  4.4  is  a  photograph  of  the  fragments  frc*n 
Tank  D.  Weights  for  each  fragment  are  listed  in  Table  4.5. 

As  may  be  seen  in  Table  4.5,  we  recovered  21  fragments 
from  Tank  D.  From  Table  4.2,  there  was  a  total  of  28  hits  on  the 
arena  wall.  Seven  of  these  fragments  remained  in  the  walls,  8  were 
found  outside  the  arena,  and  the  other  13  were  not  located.  The 
locations  where  fragments  were  found  are  shown  in  Figure  4,5. 

Tank  D  weighed  171  lbs  f 77,600  grams)  in  a 
condition  ready  for  inflation.  The  recovered  fragments  weighed  76  lbs 
(34531  grams).  Thus  the  fragments  recovered  amounted  to  45.5#  of  the 
total  tank  weight. 

4.2.5  Fragments  from  Tank  E 

A  photograph  of  the  fragments  from  Tank  E  are  shown 
in  Figure  4.6.  Table  4.6  lists  fragment  weights. 

We  recovered  25  fragments  from  Tank  E.  From  Table  4.2, 
there  were  28  hits  on  the  arena  walls.  Ten  of  these  fragments  remained 
imbedded  in  the  arena  walls,  8  were  found  outside  the  arena,  and  the 
other  10  penetrated  the  wall  at  high  velocity  and  were  not  recovered. 
The  holes  in  the  arena  walls  gave  a  general  Indication  of  the  size  of 
the  unrecovered  fragments.  Certainly  none  were  as  large  as  the  three 
largest  fragments  in  Figure  4.4  or  the  five  largest  in  Figure  4.6. 
Rather,  the  holes  gave  evidence  of  sizes  similar  to  those  represented 
by  fragments  4  through  21  for  Tank  D  and  6  through  24  for  Tank  E. 
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FIG.  4-3  FRAGMENTS  FROM  TANK  C  (2"  LINE  SPACING) 
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FIG.  4-5  APPROXIMATE  LOCATION  OF  FRAGMENTS  FOR  TANK  D 
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Locations  where  fragments  were  found  are  shown  in 
Figure  4.7.  The  numbers  inside  the  circles  correspond  to  fragment 
numbers  given  in  Table  4.6  and  Figure  4.6.  (Panel  30  was  removed  for 
tank  E. ) 

Tank  E  weighed  170  lbs  (77000  grams)  at  inflation. 

The  weight  of  the  recovered  fragments  totaled  87.7  ibs  (39,809  grams). 
Thus  the  recovered  fragments  amounted  to  51.6$  of  the  total  tank 
materials. 

4.3  Fragment  Velocities 


The  two  methods  for  measuring  fragment  velocity  are  described 
in  Section  2.3.2.  While  both  systems  did  not  operate  on  every  test, 
one  or  the  other  did  function.  Therefore,  fragment  velocities  were 
measured  on  every  test.  These  results  are  presented  below. 

4.3.1  Velocities  of  Tank  A  Fragments 


The  breakwire  system  did  not  operate  on  this  event. 
However,  the  strobe  system  yielded  adequate  pictures  of  the  fragments 
in  flight.  Times  at  which  the  strobe  flashes  occurred  are  given  in 
Table  4.7.  Figure  4.8  shows  the  strobe  picture  with  the  reflective 
backdrop.  The  shockwave  is  clearly  visible  for  the  first  three  flashes 
and  is  labeled  1-S,  2-S,  and  3-S.  The  shock  labeled  1-S  is  at  0  time, 
i.e.,  the  time  at  which  the  breakwire  system  signaled  tank  rupture. 
Obviously,  the  tank  has  burst  and  the  shockwave  formed  by  the  time  the 
wire  in  contact  with  the  tank  breaks.  At  0  +  1.5  milliseconds,  we 
photograph  the  shock  labeled  2-S  at  its  new  position.  Note  that  the 
shock  is  not  symmetrical  about  the  tank;  it  has  traveled  further  on 
the  right  where  rupture  first  occurs.  The  shock  is  not  photographed 
by  the  4th  flash  at  0  +  5  milliseconds  since  it  is  beyond  the  bounds 
of  the  scotchlite  screen  at  this  time. 

The  flashes  seen  on  Figure  4.8  and  4.9  are  self- 
illuminating  collisions  of  fragments  and  pieces  of  the  breakwire 
with  solid  objects.  Since  they  are  self-illuminating,  we  are  unable  to 
assign  a  time  to  their  occurrence.  Therefore,  we  cannot  use  them  in 
velocity  calculations. 

In  measuring  fragment  velocities,  we  assumed  that  the 
fragments  traveled  in  straight  lines.  This  may  not  be  true  since 
a  large  flat  fragment  traveling  in  a  sailing  mode  would  probably 
show  some  curvature.  However,  this  assumption  allowed  velocity 
measurements  when  the  fragment  could  not  be  identified  on  every  flash. 
All  fragments  were  assumed  to  emanate  from  the  tank's  equator 
parallel  to  the  earth’s  surface.  Fragment  locations  were  determined 
from  two  mutually  perpendicular  coordinate  systems  described  by 
camera  fields  of  view  and  the  tilt  angles  of  the  film  plane.  Tbe  tower 
camera  geometry  is  shown  in  Figure  4.10, 
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FIG.  4-7  APPROXIMATE  LOCATION  OF  FRAGMENTS  FOR  TANK  E 


FIG .  4-8  GROUND  CAMERA  PHOTO  OF  TANK  A  RUPTURE 
(SEE  SECTION  4-3-1  FOR  EXPLANATION  OF  NOTATION) 


FIG.  4-9  TOWER  CAMERA  PHOTO  OF  TANK  A  RUPTURE 
(SEE  SECTION  4-3-1  FOR  EXPLANATION  OF  NOTATION) 
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From  Figures  4.8  and  4.9  we  are  able  to  get  position- 
time  information  for  selected  fragments.  For  instance  note  the  fragment 
labeled  1-2,  1-3,  and  1-4.  The  first  digit  identifies  the  fragment. 

The  second  digit  refers  so  the  flash  number.  Thus,  1-4  refers  to 
fragment  #1  at  the  fourth  flash  or  at  5-0  milliseconds.  In  the  case 
of  the  three  fragments  chosen  for  obtaining  fragment  velocities,  all 
could  be  located  with  sufficient  accuracy  to  yield  position-time  data. 
These  measurements  are  summarized  in  Table  4.8. 

For  Tank  A,  fragments  1  and  2  were  identified  at  all 
four  strobe  flashes  in  Figure  4. '9  and  at  the  last  two  flashes  in 
Figure  4.8.  At  0  time  the  tank  is  assumed  to  be  intact.  The  apparent 
increase  in  velocity  for  fragment  2  may  be  due  to  the  fact  that  the 
fragment  perimeter  is  not  clearly  identifiable  in  the  condensation 
cloud  formed  by  the  expanding  gas.  Fragment  3  was  identifiable  on  the 
last  two  flashes  in  Figure  4.9  but  only  on  the  last  flash  in  Figure  4.8. 
Other  fragments  are  visible  in  Figure  4.9  tut  cannot  be  found  on 
Figure  4.8. 

4.3.2  Velocities  of  Tank  B  Fragments 

Both  the  strobe  system  and  the  breakwire  system 
yielded  velocity  data  for  this  event.  Figure  4.11  shows  that  tank 
rupture  has  occurred  and  a  shock  marked  1-S  is  emerging  from  the 
upper  end  of  the  tank  when  the  start  signal  occurs.  (The  shock  is 
also  seen  on  the  second  and  third  flashes  and  is  marked  2-S  and  3-S 
respectively.)  At  this  time  the  distance  between  the  tank  perimeter 
and  the  stop  wire  support  is  0.77  ft.  The  counter-measured  time 
between  the  breakwire  signals  was  0.632  milliseconds  giving  an  average 
fragment  velocity  of  1215  ft/sec.  This  value  is  not  out  of  line  with 
velocities  measured  by  the  photosystem. 

Fragment  velocity  measurements  were  made  from  Figures 
A. 11  and  4.12.  Unfortunately  the  2  fragments  labeled  2-4,  2-3,  1-4,  and 
1-3  in  Figure  4.12  could  only  be  identified  for  the  third  flash  in 
Figure  4.11.  From  this  we  were  able  to  obtain  a  trajectory  and  there¬ 
fore  some  velocity  information. 

Fragment  1  travels  a  distance  of  4.56  ft  in  4.5  milli¬ 
seconds  for  an  overall  velocity  of  1010  ft/sec.  Its  travel  time 
between  the  last  two  flashes  is  1.5  milliseconds.  It  covered  the 
distance  between  3-06  and  4.56  ft  or  1.5  ft.  Average  velocity  over  this 
interval  is  1000  ft/se^. 

Fragment  2  traveled  a  distance  of  4.65  ft  over  4.5 
milliseconds  an  overall  velocity  of  1030  ft/sec.  Its  travel  time 
between  the  last  two  flashes  was  1.5  milliseconds.  During  this  time 
it  covered  a  distance  of  1.35  ft  for  a  velocity  of  900  ft/sec.  These 
measurements  are  summarized  in  Table  4.8. 
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FIG.  4-12  TOWER  CAMERA  PHOTO  OF  TANK  B  RUPTURE 
(SEE  SECTION  4-3-1  FOR  EXPLANATION  OF  NOTATION) 
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4.3,3  Velocities  of  Tank  C  Fragments 


Fragment  velocities  for  this  tank  were  obtained  from 
both  the  breakwire  system  and  the  strobe  system.  However,  data  from 
both  systems  require  interpretation.  For  the  breakwire  system,  the 
two  electronic  counters  gave  different  readings.  One  counter  gave 
a  travel  time  for  12.5  inches  of  O.876  milliseconds  for  a  velocity 
equal  to  1190  ft/sec.  The  other  gave  a  travel  time  of  0.773  milli¬ 
seconds  over  the  same  distance  for  a  velocity  of  1350  ft/sec. 

The  shutter  for  the  ground  camera  did  not  close 
before  the  photofloods  came  on  to  illuminate  the  dummies.  Therefore 
the  film  from  this  camera  was  overexposed  (Fig.  4.13).  While  the 
shadows  left  by  the  shockwave  were  not  washed  out,  we  can  identify 
only  one  fragment.  Therefore,  a  fragment  path  is  based  on  this  and 
Tig.  4.14  to  obtain  velocity.  From  fragment  1  we  get  a  total  distance  tra¬ 
veled  of  5.20  ft  in  4.50  milliseconds  or  an  average  velocity  of  1150  ft/sec. 
Over  the  first  3.00  milliseconds,  it  traveled  3.60  ft  or  1200  ft/sec. 

Over  the  final  interval  it  traveled  1020  ft/sec. 

4.3.4  Velocities  of  Tanks  D&E  Fragments 

Both  the  strobe  system  and  the  breakwire  system 
functioned  during  the  rupture  of  these  two  large  vessels.  Because 
of  the  condensation  cloud,  no  fragments  could  be  seen  in  the  photo¬ 
graphs.  Therefore,  only  the  breakwire  system  yielded  fragment 
velocities. 


For  Tank  D  the  distance  traveled  was  0.97  -  0.062  ft 
over  the  measured  time  interval  of  0.689  milliseconds.  The  0.062  ft 
is  the  accuracy  we  could  maintain  of  the  cage  distance  from  the  tank. 
This  gives  an  average  velocity  of  1400  ft/sec  over  this  interval. 

The  reason  that  no  fragments  could  be  clearly  identi¬ 
fied  in  the  strobe  photographs  was  that  condensation  at  the  front  of 
the  fill  gas  out  ran  and  obscured  everything  behind  it. 

On  Tank  D,  the  strobe  was  flashed  at  0,  0  +  1.5, 

0  +  3.0,  and  0  +  4.5  milliseconds.  Figure  4,15  shows  that  the  cloud 
has  reached  a  distance  of  about  7  ft  (distance  to  the  first  gages) 
in  4.5  milliseconds. 

Figure  4.l6  shows  the  shock  at  0  +  1.5  milliseconds. 

No  fragments  are  visible.  However,  the  shock  front  has  covered  a 
distance  of  about  3  ft  in  the  vertical  direction  during  this  period, 
an  average  velocity  of  2,000  ft/sec.  (Distances  are  given  from  the 
tank  perimeter  to  the  shock  front.)  To  che  right  in  Figure  4.16,  the 
shock  has  traveled  less  than  2  ft  for  a  velocity  of  1300  ft/sec. 

No  fragments  are  visible  ahead  of  the  condensation  front.  The  multiple 
shock  structure  at  the  right  may  be  due  to  fragment  bow  waves  but  is 
more  probably  due  to  the  jetting  of  the  N  at  tank  rupture. 

9 


78 


FIG.  4-13  GROUND  CAMERA  PHOTO  OF  TANK  C  RUPTURE 
{SEE  SECTION  4-3-1  FOR  EXPLANATION  OF  NOTATION) 
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FIG.  4-14  TOWER  CAMERA  PHOTO  OF  TANK  C  RUPTURE 
(SEE  SECTION  4-3-1  FOR  EXPLANATION  OF  NOTATION) 
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The  strobe  system  was  changed  for  Tank  E  in  an 
attempt  to  photograph  fragments.  Height  of  the  scotchlite  screen 
for  the  ground  camera  was  increased  to  12  ft  and  width  was  increased 
to  24  ft.  No  scotchlite  was  used  in  the  center  of  the  screen  since 
it  was  expected  that  the  condensation  cloud  would  cover  this  area. 

Hie  strobe  times  were  changed  to  flash  at  0,  0  +  6,  0  +  7,  and  0  +  8 
milliseconds.  It  was  hoped  that  the  fragments  would  out-distance 
the  cloud  during  this  extended  period.  Fragments  could  be  clearly 
identified  above  the  cloud  in  Figure  4.17,  the  tower  camera  photograph 
but  not  in  Figure  4.18,  the  ground  camera  photograph.  Therefore, 
trajectory  information  was  not  available  and  no  photo  velocities  could 
.be  calculated.  Considerable  flashing  is  evident  in  Figure  4.17 
probably  due  to  fragments  encountering  bits  of  the  wire  cage  and  other 
solid  objects. 
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FIG 4-18  GROUND  CAMERA  PHOTO  OF  TANK  E  RUPTURE 
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5.  Anthropomorphic  Dummies 


For  purposes  of  producing  training  films,  anthropomorphic  dummies 
were  placed  near  the  rupturing  tanks.  Motion  picture  photography 
recorded  the  visual  effects  of  the  "blast  waves  on  the  dummies. 

For  tank3  A,  B,  C,  and  D  the  three  dummies  were  placed  as  shown  in 
Figure  2.1.  Tank  to  dummy  distances  for  these  four  tests  were  8.5, 
12.8,  and  13.7  feet.  For  the  tank  E  test,  the  dummies  were  placed 
at  slightly  different  relative  positions.  For  this  event,  tank  to 
dummy  distances  were  10.0,  l4.5>  and  15.2  feet.  The  dummies  were 
supported  in  an  upright  position  by  1”  diameter  aluminum  pipes  driven 
into  the  ground  so  that  their  tops  were  at  waist  level.  They  were 
placed  in  a  manner  such  that  the  pipes  offered  little  support  against 
the  blast. 

None  of  the  dummies  were  knooked  over  by  the  blast  from  any  tank 
rupture,  though  some  were  knocked  over  by  the  collapse  of  the  arena 
walls.  For  events  A,  B,  and  C,  the  closest  dummy  was  in  an  area  where 
a  man  would  suffer  ear  drum  damage  —  side-on  overpressures  above  5.0 
psi.  For  tanks  D  and  E,  all  the  dummies  were  in  an  area  where  ear 
drum  damage  would  definitely  occur.  The  closest  dummy  on  event  D  was 
in  the  25-psi  region,  a  level  where  the  onset  of  lung  damage  occurs. 

It  is  unlikely  that  death  would  occur  from  blast  alone.  The  frag¬ 
ments,  however,  could  have  proven  lethal  if  a  hit  occurred  in  the 
right  spot.  On  tank  A,  the  closest  dummy  was  hit  in  the  buttocks  by 
a  large  fragment,  tearing  out  the  seat  of  his  coveralls.  In  addition 
to  embarrassment,  a  man  would  have  suffered  severe  lacerations.  Another 
fragment  from  tank  A  sliced  into  and  remained  in  the  leg  of  the  second 
dummy.  A  man’s  leg  would  have  been  severely  cut,  with  a  probable  bone 
fracture.  A  fragment  from  tank  D  struck  the  leg  of  the  second  dummy 
edge  on  with  sufficient  force  to  bend  the  heavy  metal  I-beam  in  the 
leg  In  a  right  angle.  A  man's  leg  would  have  been  severed. 

All  the  tanks  threw  up  dirt  and  gravel  with  sufficient  force  to 
penetrate  a  man's  skin  and  to  cause  severe  eye  damage  and  probable 
blindness . 
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6.  Discussion 

6.1  Tank  Pressurization 

Tank  pressurization  results  showed  that  none  of  the  vessels 
burst  at  less  than  their  rated  burst  pressure.  Hie  fuel  tanks  rated 
at  460  psig  burst  pressure  actually  burst  at  625  psig  for  the  aerozine 
fuel  tank  and  at  600  psig  for  the  oxidizer  tank.  Hie  small  helium 
tank,  rated  at  7>500  psig  actually  ruptured  at  8,000  psig.  The  two 
large  helium  vessels  were  rated  at  8,000  psig  burst  pressure.  The 
first  of  these  vessels  burst  at  exactly  8,000  psig.  The  second  tank 
was  pressure  cycled  55  times.  This  tank  burst  at  8,130  psig. 

The  tanks  were  radiographed  prior  to  the  test3  to  locate 
voids,  flaws,  sub-standard  welds,  or  other  defects.  Dyepenetrate 
checks  were  also  performed.  No  deficiencies  were  found.  Indeed, 
tanks  ruptured  in  what  seemed  a  random  manner.  Hiere  were  no  patterns 
such  as  rupture  along  weldments. 

Hiere  was  no  apparent  residual  effect  from  cycling  the  larg^ 
helium  tank.  It3  burst  pressure  wa3  slightly  but  not  significantly 
higher  than  that  of  the  similar  but  uncycled  tank.  Hie  pattern  of 
fragmentation  seemed  unaffected.  About  the  same  total  weight  and 
number  cf  fragments  were  recovered  from  both  the  large  helium  tanks 

6.2  Blast  Measurements 

Side-on  and  face-on  pressure  measurements  showed  a  lack  of 
symmetry  about  GZ  for  the  shock  wave  generated  by  tank  rupt*a*e.  Hie 
slope  of  the  pressure  distance  curve  was  less  for  tank  rupture  data 
than  for  TNT  data.  For  this  reason,  a  single  TNT  equivalent  does  not 
characterize  the  TNT  tank  rupture  data.  Even  so,  tank  ru yture  blast 
is  a  function  of  rupture  pressure  and  tank  volume.  Therefore,  we 
expect  tank  rupture  airblast  to  scale  geometrically,  fhe 
pressure  distance  curves  from  such  tanks  coincide  if  distances  are 
given  3n  terms  of  tank  radii.  This  is  cube  root  scaling. 

Shock  overpressures  fran  tanks  C,  D,  and  E  have  been  averaged 
and  are  plotted  as  a  function  of  tank  radii  in  r igure  6.1.  Therefore, 
this  curve  generalizes  the  pressure  distance  curve  for  similar  tanks 
of  any  size  burst  at  8000  psig.  Caution  is  recommended  since  tank  rup¬ 
ture  5s  asymmetric  and  shock  overpressures  tnat  are  30 $  higher  than 
those  in  Figure  6.1  may  be  encountered. 

The  TNT  curve  is  based  on  the  TNT  weight  whose  blast  energy 
is  equal  to  the  tank  rupture  energy  fo.-  the  tank  size  chosen  where  the 
fill  gas  is  air  or  Na .  The  TNT  energy  is  taken  as  1018  calories/gram. 
Based  on  the  two  curves  of  Figure  6  1,  the  equivalent  weight  of  tank 
rupture  in  terms  of  TNT  has  been  calculated  and  is  plotted  as  a 
function  of  tank  radius  in  Figure  6.2. 
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FIG.  6.1  PEAK  AIRBLAST  OVERPRESSURES  FROM  AN  8,000  PSIG 
TANK  RUPTURE  VS  DISTANCE  IN  TANK  RADII 
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FIG.  6.2  TNT  EQUIVALENCE  OF  THE  AIR  BLAST  FROM  AN  8,000  PSIG 
TANK  RUPTURE  BASED  ON  PEAK  SHOCK  OVERPRESSURE 
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The  positive  shockwave  impulse  from  tanks  C,  D,  and  E  have 
been  combined  and  are  shown  in  Figure  6.3.  Hie  curve  is  universal 
for  similar  spherical  tanks  burst  at  3, GCO  psig.  Again  deviations 
Tron  the  curve  may  be  50$  due  to  tank  rupture  blast  wave  asymmetry. 

,:e  assumed  that  tan!:  rupture  equivalent  'weight  based  on  positive  shock 
wave  impulse  was  constant  based  on  TNT.  From  Figure  6.3*  E V/  for  tank 
rupture  is  about  C.Jo  x  EW  for  TNT  when  the  T1W!  weight  is  calculated 
on' the  basis  of  isentropic  expansion. 

The  comparison  of  the  airblast  data  from  tanks  C,  D,  and  E 
with  that  from  the  TNT  charges  is  correct  in  that  both  the  tanks  and 
TNT  charges  were  spherical  in  shape.  The  comparison  of  the  airblast 
data  from  tanks  A  and  B  with  the  spherical  TNT  charge  data  may  be 
questioned  since  these  tanks  were  cylindrical.  Certainly  the  blast 
field  generated  by  a  spherical  TNT  charge  would  differ  from  that 
generated  by  a  TNT  cylinder  of  equal  weight.  (Reference  (8))  Moreover, 
the  pressure  field  generated  by  a  cylindrical  charge  is  dependent  on 
the  length  to  diameter  ratio  (L/D)  and  the  point  of  detonation.  Kie 
point  of  detonation  determines  the  pattern  of  the  detonation  wave 
propagation  which  in  turn  determines  the  blast  field  about  the  cylinder. 

In  the  case  of  the  cylindrical  tanks,  there  is  no  detonation 
wave.  Therefore,  comparing  these  airblast  data  with  that  from  a  cylin¬ 
drical  TNT  charge  has  no  more  validity  than  comparing  them  with  spheri¬ 
cal  TNT  data.  The  similarity  of  the  relationship  between  the  cylindri¬ 
cal  tank- spherical  charge  and  the  spherical  tank- spherical  charge  data 
is  seen  in  Figures  3.8  through  3.24. 

Even  if  we  used  cylindrical  charges  of  the  same  scaled  shape, 
the  blast  field  would  still  depend  on  the  detonation  scheme.  Since 
we  had  no  physical  reason  for  choosing  any  one  detonation  scheme  over 
another,  no  cylinders  were  fired. 

Since  the  L/D  of  the  cylindrical  tanks  was  small  (about  2/1) 
and  the  tanks  had  hemispherical  ends,  we  treated  them  as  equivalent 
volume  spheres .  To  generalize,  we  plotted  the  mean  peak  shock  over- 
pressurer  vs  distance  in  tank  radii  for  the  two  cylindrical  l 

tanks  in  Figure  6.4.  (The  radius  is  that  of  a  sphere  v?hose  volume  is 
equal  to  that  of  a  cylinder.)  The  pressure  distance  curve  in 

Figure  6.4  should  be  representative  of  that  from  similar  tanks  burst 
at  600  psig  -  with  allowances  of  about  3C$  for  blast  wave  asymmetry . 

Based  on  the  curves  in  Figure  5.4,  tank  rupture  blast  equiva¬ 
lence  based  on  peak  shock  overpressure  has  been  calculated  and  is 
shown  as  a  function  of  tank  radiis  in  Figure  6.5. 

'Hie  positive  shockwave  impulse  from  tne  500  psig  tank  ruptures 
have  been  combined  and  are  shown  as  a  function  of  tank  radius  in 
Figure  6.6.  Again,  this  curve  should  be  representative  of  the  positive 
impulse  from  any  sized  similar  tanks  that  burst  at  600  psig.  From 
Figure  6.6,  the  tank  rupture  TNT  equivalence  in  terms  of  the  weight 
of  TNT  as  determined  from  isentropic  expansion  energy  is  considered  a 
constant.  In  this  case,  the  tank  rupture  EM/TNT  EVJ  averages  about 
1.73* 
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FIG.  6.3  POSITIVE  SHOCKWAVE  IMPULSE  FROM  AN  8,000  PSIG 
TANK  RUPTURE  VS  DISTANCE  IN  TANK  RADII 
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FIG.  6.4  PEAK  AIRBLAST  OVERPRESSURE  FROM  A  600  PSIG  TANK  RUPTURE 
VS  DISTANCE  IN  TANK  RADII 
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FIG.  6.5  TNT  EQUIVALENCE  OF  THE  AIRBLAST  FROM  A  600  PSIG 
TANK  RUPTURE  BASED  ON  PEAK  SHOCK  OVERPRESSURE 
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FIG.  6.6  POSITIVE  SHOCK  WAVE  IMPULSE  FROM  A  600  PSIG 
TANK  RUPTURE  VS  DISTANCE  IN  TANK  RADII 
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Using  rhe  results  of  these  experiments  it  is  possible  to 
predict  the  blast  parameters  generated  by  the  rupture  of  similar  tanks 
of  different  volumes  providing  rupture  pressures  are  equal.  It  is  not 
simple  to  predict  the  blast  field  for  different  rupture  pressures. 
However,  slight  modifications  to  computer  programs  now  being  used  at 
NOL  would  allow  such  predictions.  The  experimental  data  from  this 
program  would  serve  as  a  check  on  the  validity  of  such  calculations . 

6 . 3  Fragmentation 

6.3.1  Fragment  Density 

Procedures  used  for  obtaining  fragment  size  and  distribution 
from  munitions  is  misleading  in  the  case  of  the  tanks.  The  number  of 
fragments  from  even  a  small  bomb  are  orders  of  magnitude  greater  than 
the  number  from  the  tanks .  Statistical  treatments  for  munitions 
fragments  are  justifiable.  In  the  case  of  the  tanks,  there  were  only 
three  hits  on  a  circular  arena  from  one  tank:  this  hardly  leads  to  a 
statistical  analysis . 

A  further  complication  in  the  tank  case  was  the  ricocheting 
fragments .  In  munitions  work,  a  barricade  is  erected  to  intercept 
such  fragments .  In  our  case,  we  were  unable  to  take  this  precaution 
because  a  barricade  would  have  ruined  blast  measurements .  A  better 
technique  is  to  estimate  the  total  number  of  fragments  based  on  the 
number  and  weight  of  those  recovered. 

In  the  case  of  tanks  A  and  B,  the  number  of  heavy  tank  parts 
recovered  could  bias  this  estimate  since  with  their  weight  added,  we 
recovered  about  70$  of  the  tanks  by  weight.  Since  all  the  heavy 
parts  were  recovered,  the  unrecovered  fragments  consisted  of  the  tank 
walls  only. 

We  actually  recovered  56$  of  the  material  making  up  the 
tank  walls.  Of  the  36  fragments  recovered,  32  were  of  tank  wall  ma¬ 
terial.  On  this  basis  32  fragments  equal  56$  of  the  material,  there¬ 
fore  the  total  number  of  fragments  estimated  from  tank  A  would  be  57- 
With  the  heavy  fragments,  this  makes  a  total  of  6l.  We  further  assume 
that  these  fragments  were  scattered  equally  over  a  360°  solid  angle. 

For  tank  B,  6l$  of  the  tank  wall  material  was  recovered. 

'Bils  was  divided  into  33  fragments.  The  estimated  total  number  of 
tank  wail  fragments  then  would  be  5^.  Including  the  heavj  tank  parts, 
the  total  is  57. 

There  is  a  possibility  that  in  some  cases,  the  heavy  wire 
frame  supporting  the  second  breakwire  may  have  cut  some  of  the  larger 
fragments  into  many  smaller  pieces.  Thus  had  the  tank  ruptured  in  an 
unencumbered  environment,  fewer  fragments  would  have  resulted. 

For  tank  C,  because  of  the  richcchet  problem,  only  21$  of 
the  tank  was  recovered.  This  consisted  of  8  fragments.  A  wild  esti¬ 
mate  of  the  total  number  of  fragnents  from  tank  C  would  be  38  fragments. 
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For  tank  D  we  recovered  21  fragments,  45$  of  the  tank  by 
weight.  Based  on  this,  47  fragments  were  produced.  For  tank  E  we 
recovered  25  fragments  making  up  51.6$  of  the  tank’s  weight.  Based 
on  this  the  total  number  of  fragments  from  tank  E  would  have  been  48. 

In  the  above,  the  flange  (fragment  25  in  Figure  4.6)  has  been  ignored 
since  :.ts  accounting  would  not  influence  the  estimate. 

Of  the  21  fragments  recovered  from  Tank  u,  six  were 
inside  the  arena.  The  other  15  were  either  in  the  arena  wall  or  just 
outside  the  walls.  From  Table  4.2  there  were  28  arena  wall  hits,  Hie 
above  accounts  for  15  hits .  Thus  there  were  13  more  fragments  with 
trajectories,  either  direct  or  by  ricochet*  that  intercepted  the  arena 
wall.  Assume  that  all  the  fragments  that  were. recovered  had  tra¬ 
jectories  that  were  below  the  arena  walls.  (This  means  that  no 
fragments  had  trajectories  near  enough  to  the  vertical  to  fall  back  in 
the  arena.)  With  this  assumption,  we  add  the  13  fragments  not  found 
to  those  21  found  for  a  total  of  34  fragments .  This  means  that  the 
trajectories  of  all  34,  or  72$  were  below  the  arena  walls .  Yet  the 
top  of  the  arena  walls  defined  a.i  angle  that  covered  5S$  of  the  tank's 
surface . 

Based  on  a  similar  argument  for  tank  E,  35  fragments  or  73 $ 
had  trajectories  below  the  top  of  the  arena  walls .  Again  the  angle 
described  by  the  arena  walls  intercepted  only  59$  of  the  tank’s 
surface . 

We  are  left  with  these  conclusions:  (a)  Our  method  of 
estimating  number  of  fragments  is  incorrect,  (b)  The  fragment  dis¬ 
tribution  is  not  homogeneous  over  a  solid  angle  described  by  a  spherical 
shell,  (c)  All  the  fragments  found  on  the  floor  of  the  arena  went 
straight  up  and  fell  back  in  the  arena.  (Fragments  15*  17*  18,  19* 

22,  23,  and  25  in  Figure  4.7  for  instance.)  1110  most  probable  factor 
is  that  we  are  underestimating  the  number  of  fragments  because  of  the 
influence  of  the  fragment  numbered  one  in  Tables  4.5  and  4.6  on 
fragment  distribution  and  average  weight. 


To  aid  in  experimental  planning,  a  short  computer  program 
was  devised  to  give  "ball  park"  estimates  of  fragment  velocities . 

Tliese  calculations  are  presented  In  Appendix  B.  Calculated  fragment 
velocities  tended  to  be  20  to  40$  lower-  than  those  measured  by  the 
breakwire  and  photographic  system. 

Ihe  calculated  values  for  the  cylindrical  tanks  (tanks  A  and 
B)  in  Figure  5  of  Appendix  B  were  even  lower.  However,  a  later  calcu¬ 
lation  using  the  correct  burst  pressure  and  density  for  the  tank  wall 
material  yielded  maximum  fragment  velocity  closer  to  the  measured 
value.  (Figure  6.7)  The  rapid  decay  of  the  calculated  velocity  with 
time  results  from  the  model  used  for  the  calculations ,  In  the  model, 
the  fragments  are  tied  to  the  gas  front  in  an  arrangement  somewhat  as 
though  they  were  on  a  rubber  band.  The  significant  thing,  from  the 
standpoint  of  damage,  is  that  the  fragments  do  not  slow  down  rapidly. 
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In  Figures  4.8  and  4.11,  tank  rupture  is  seen  to  have  already- 
taken  place  by  the  time  the  first  flash  occurs .  The  counter  for 
velocity  measurements  started  at  the  same  time  that  the  first  flash 
occurred,  within  5  microseconds.  The  question  then  arises,  were  frag¬ 
ments  already  in  flight  when  the  counters  started?  If  so,  then  the 
breakwire  system  would  yield  a  velocity  greater  than  actually  existed. 

For  tank  A,  the  breakwire  system  did  not  operate.  Therefore, 
we  are  unable  to  make  a  comparison.  Note  however,  in  Table  4.8,  the 
average  fragment  velocity  over  the  first  1.5  milliseconds  of  time  as 
measured  by  the  strobe  system  was  slightly  over  1000  ft/sec. 

For  Tank  B  the  breakwire  velocity  was  1215  ft/sec.  The 
strobe  system  did  not  yield  fragment  velocities  over  the  first 
1.5  millisecond  time  interval.  However,  velocities  measured  at  later 
times  by  the  strobe  system  were  not  greatly  different  from  those 
of  Tank  A.  For  Tank  C,  the  1200  ft/sec  fragment  velocity  measured 
by  the  strobe  system  over  the  first  3.0  milliseconds  is  not  out  of 
line  with  the  1270  +  80  ft/sec  measured  by  the  breakwire  system. 

For  Tanks  D  and  E,  the  strobe  system  did  not  yield  fragment 
velocities.  However  from  Figure  4.15  we  can  put  an  upper  velocity 
limit  on  the  Tank  D  fragments.  In  Figure  4.15,  no  fragments  are 
visible  outside  the  condensation  cloud  at  0  +  4.5  milliseconds.  Since 
the  cloud  is  at  a  7  ft  radius  at  this  time  and  the  tank  radius  is 
1.13  ft,  the  distance  for  a  fragment  to  travel  would  be  5.87  ft. 
Therefore,  no  fragment  has  traveled  at  an  average  velocity  greater 
than  1300  ft/sec  over  the  first  4.5  milliseconds.  In  Figure  4.17, 
fragments  from  Tank  E  are  clearly  outside  the  condensation  cloud  at 
0  +  S  milliseconds.  Since  the  cloud  is  about  30.5  ft  in  diameter  at 
this  time,  the  fragments  must  have  traveled  a  distance  greater  than 
9.4  ft  from  the  tank  perimeter.  Therefore  the  fragment  velocity  was 
greater  than  1170  ft/sec  averaged  over  8  milliseconds . 

6 . 4  Qualitative  Damage 


It  was  shown  in  Section  5  that  the  fragment  hazard  to  person¬ 
nel  is  severe  from  any  of  the  tanks .  Certainly  people  deployed  as 
were  our  anthropomorphic  dummies  would  have  been  maimed,  possibly 
killed.  The  damage  potential  of  heavy  fragments  —  several  pounds  — 
traveling  at  supersonic  velocities  is  awesome.  These  fragments  sheared 
off  heavy  3”  pipes,  passed  through  6!1  of  panel  broadside  on,  and  were 
found  over  300  yards  from  Ground  Zero! 

Even  the  blast  from  any  of  the  tanks  would  rupture  eardrums 
and,  in  the  case  of  tanks  D  and  E,  cause  lung  damage.  More  impressive, 
the  blast  from  tanks  D  and  E  seemed  to  do  things  that  would  not  be 
expected  from  10  pounds  of  TNT.  Specifically  the  blast  from  tank  D 
drove  the  4*  x  4*  x  1.5’  concrete  firing  pad  about  4  inches  into  the 
ground.  The  tank  E  rupture  drove  the  pad  another  5  inches  into  the 
ground.  Both  tanks  caved  in  and  collapsed  the  arena  walls.  The  8.2 
pounds  of  TNT  neither  caved  in  the  wall  nor  displaced  the  concrete 
pad  any  further,  nor  would  the  TNT  be  expected  to  do  so. 
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Side- on  overpressures  and  positive  shock  wave  durations  from 
tanks  D  and  E  were  different  from  those  for  a  TNT  equivalent  charge 
based  on  isentropic  expansion;  but  they  combined  to  give  positive 
shock  wave  impulses  quite  comparable  to  those  for  the  TNT.  Neverthe¬ 
less,  blast  from  tank  rupture,  qualitatively  at  least,  seemed  more 
violent  than  that  from  TNT  at  comparable  distances. 
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7 .  Conclusions 

Our  studios  proved  that  no  single  TNT  equivalent  can  be  used  to 
describe  the  blast  generated  by  a  rupturing  pressure  vessel.  For 
instance,  side-on  pressures  from  tank  rupture  are  lower  than  those 
for  an  equivalent  TNT  charge  at  blast  pressures  above  10  to  15  psi 
but  are  higher  than  those  for  TNT  below  10  to  15  psi.  However,  the 
tank-rupture  blast  pressures  combined  with  the  positive  shock  wave 
durations  yield  positive  shock  wave  impulse  values  whose  impulse- 
distance  relation  was  similar  in  slope  to  that  for  TNT.  In  the  case 
of  the  large,  high  pressure  vessels,  impulse  for  tank  rupture  and  for 
the  TNT  equivalent  were  not  significantly  different  quantitatively. 

The  large,  high  velocity  fragments  generated  by  tank  rupture  will 
prove  to  be  a  significant  hazard.  Tanks  A  and  B,  similar  cylindrical 
tanks,  broke  into  about  the  same  number  of  fragments  and  generated 
similar  fragment  velocities.  The  same  reproducibility  of  number  of 
fragments  and  similar  velocities  were  recorded  for  the  high  pressure 
spherical  vessels,  tanks  D  and  E. 

We  can  make  a  general  comparison  of  blast  and  fragment  parameters 
generated  by  tank  rupture  and  an  equivalent  TNT  charge.  Side-on  pres¬ 
sures  are  higher  for  TNT  above  10  to  6  psi  and  lower  for  TNT  at 
pressures  below  these  values.  Face-on  tank  pressures  show  a  similar 
relati  n  to  face-on  TNT  pressures.  Positive  shock  wave  durations 
are  Ion"--  *r  for  tank  rupture  than  for  TNT.  Impulse,  both  face-on 
and  side-on,  is  similar  for  TNT  and  tank  rupture.  Damage,  depending 
on  distance,  may  be  greater  for  tank  rupture.  Tank-rupture  fragments 
were  larger  than  would  be  expected  from  a  cased  TNT  charge.  Fragment 
velocities  would  be  higher  for  a  cased  TNT  charge  than  for  tank 
rupture . 

Further  analysis  of  the  tank  rupture  data  and  phenomena  are 
required  to  correlate  damage  with  the  blast  parameters. 

The  work  reported  here  can  be  extended  by  calculations  to  other 
tanks  with  different  shapes,  volumes,  and  rupture  pressures.  The 
computational  efforts  in  the  appendices  can  be  refined  to  handle  such 
work . 
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Table  3.1 

Blast  data  from  24"  Cylinder;  Tank  A 
Burst  Pressure  =  625  psi 
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Gage  to  tank  center  distance  in  feet 
Maximum  shockwave  pressure  in  psig 
Arrival  time  in  milliseconds 

Positive  shockwave  duration  in  milliseconds  (ms) 
Positive  shockwave  Impulse  in  psi-ms 
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Table  4.1 


Fragmentation  Data  from  Tank  A 


Fragment 

Weight  Where 

Fragment 

Weight 

Where 

Number 

Grams  Found 

Number 

Grams 

Found 

( Panel ) 

(Panel) 

1 

1248  (note  1) 

22 

6.4 

2 

177  6 

23 

5.5 

3 

149  # 

24 

5.2 

4 

109  (note  2) 

25 

5.2 

5 

105  2 

2b 

5.0 

6 

100 

27 

3.7 

7 

96 

28 

3.4 

8 

78  (note  2) 

29 

2.7 

12 

9 

53  5 

30 

2.5 

10 

53 

31 

2.2 

13 

11 

47 

32 

2.0 

12 

46 

33 

1.9 

13 

42 

34 

145. 

14 

39 

35 

89. 

15 

33 

36 

102. 

16 

18  4 

Total 

55377“ 

17 

15 

18 

14  1 

IS 

13  4 

20 

11 

21 

8.7(note  3) 

Note  1: 

Found  at  GZ 

Note  2: 

Found  in  camera  shield  at  arena 

opening 

Note  3: 

Found  in  leg  of  dummy 

*  Total  tank  weignt  =  3850  grams 
’  Tank  wall  thickness  =  0.018''  to  0.024" 


*  Material  density  =  4.4l  gm/cc 
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Table  4.2 
Panel  Hits 


mANK 


Panel 

A 

B 

C 

D 

E 

1 

1 

1 

2 

2 

2 

1 

1 

-j 

2 

4 

2 

1 

R 

1 

6 

•» 

j. 

■} 

-J 

7 

g 

3 

Q 

1 

2 

1C 

1 

11 

12 

2 

2 

13 

2 

1 

1 

•3 

14 

15 

1 

1 

.1 

16 

2 

1 

17 

1 

18 

4 

19 

1 

20 

1 

21 

1 

22 

2 

2 

23 

p 

24 

1 

3 

25 

2 

1 

26 

1 

2 

1 

27 

1 

1 

28 

1 

29 

2 

30 

1 

3 

31* 

? 

Total 

12 

15 

n 

0 

rjO 

28 

Recovered 

36 

36 

Q 

o 

21 

25 

*  31  is  prote 

ctive 

;ood  shield 

in 

front 

of  came 

openj  ng 


arena 
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Fragment 

Number 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22A 
22 
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Table  4.3 


Fragmentation  Data  from  Tank  B 


Weight 

Where 

Fragment  Weight 

Where 

Grams 

Found 

Number 

Grams 

Found 

( Panel ) 

(Panel) 

1503  (Note  1) 

23 

10 

397 

24 

6.6 

142 

1 

25 

6.4 

125 

26 

6.2 

114 

30 

27 

4.6 

100 

24 

28 

3.2 

92 

23 

29 

2.7 

72 

30 

2.0 

65 

31 

1.1 

63 

5 

32 

0.4 

48 

33 

0.3 

44 

34 

101. 

38 

35 

70. 

33 

25 

Total 

3300 

33 

32 

31 

31 

oft 

Note  1: 

Found  at  GZ 

26 

f. 

*  Total 

tank  weight  = 

4100  gram 

oh 

C 

Tank  wall  thickness 

=  0.018" 

*-*T 

23 

0 

25 

*  Material  density  =  4.4l  gm/cc 

20 

6 

.024" 
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Fragment 
Numbe  r 

1 

2 

3 

4 

5 


7 

8 


Table  4.4 

Fragmentation  Data  .from  Tank 

Weight  Where 

Found 

(Grama)  (Panel) 

300  (Nbte  1) 

167  20 

67  (Note  2) 

30.9 

12.9 

13.1  (Note  3) 

15.6 

12.6 


Total 


617.1 


Note  1:  Found  at  GZ 

Note  2:  Through  scotchlite  at  Panel  13 
Note  3:  Through  scotchlite  at  Panel  15 


Total  tank  weight  =  2860  grams 
Tank  wall  thickness  =  0.108" 
Material  density  =  4.44  gm/cc 
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Table  4.5 


Fragmentation  Data  from  Tank  D 


Fragment  Weight  where 
Number  (Grams)  Found 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


7321 

2923 

2639 

2497 

2270 

2156 

2156 


2157 

2043 

2015 

1731 

3731 


482 


454 

227 

199 

85 

57 

111 


GZ 

Outside  arena  panel  3 
Outside  arena  panel  9 
Outside  approximately  300  yards 
4'  from  GZ  tov;ard  panel  27 

panel  10  low  outside 


Panel  25  outside 
Panel  18  outside 
256  yards  outside  over  panel  26 

Panel  18 
Panel  30  low 
Panel  10 
Base  panel  2 
Panel  19 
Panel  12 
Panel  18 

10'  from  GZ  toward  panel  19 
17 '  from  GZ  toward  panel  7 
14 '  from  GZ 
Base  panel  2 
Panel  27 


Total  weight  34,531  gms  =  76  lbs 
Average  thickness  0.362  inches 
Density  4,4 6  g/cc 
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Table  4.6 


Fragmentation  Data  from  Tank  E 


Fragment 

'/eight 

Where 

Number 

(Grams ) 

Found 

1 

8143 

Panel  18  high 

2 

3802 

Panel  3 

3 

3036 

Outside  arena  behind  panel  1 

4 

2610 

Outside  arena  near  camera  station  2 

5 

2412 

Panel  8  found  at  base 

6 

2270 

Panel  21  high 

7 

2270 

Panel  20  low 

8 

1873 

Panel  17  low  found  outside 

9 

1759 

Panel  26  low  found  outside 

10 

1674 

Panel  15 

11 

1617 

Panel  19  high  found  outside 

12 

1277 

Panel  3 

13 

1022 

Panel  27  middle 

14 

1050 

Panel  25  low 

15 

738 

Inside  arena 

16 

624 

Panel  8 

17 

539 

Inside  arena 

18 

454 

Panel  27  "base  on  ground 

19 

397 

Inside  arena 

20 

199 

Panel  8 

21 

170 

Panel  18 

22 

142 

Inside  arena 

23 

61 

Inside  arena 

24 

38 

Panel  29  high 

25 

1532 

GZ  intake  block 

Total  weight 

39,809 

gms  =  87.7  lbs 

Average  thickness  0.365  inches 

Density  4.46  g/cc 
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Table  4.7 

Times  of  Strobe  Plashes 

(All  times  given  in  milliseconds  after  tank  rupture) 


Tank 

1st  Flash 

2nd  Flash 

3rd  Flash 

4th  Flash 

A 

0 

1.50 

3.00 

5.00 

B 

0 

1.50 

3.00 

4.50 

C 

0 

1.50 

3.00 

4.50 

D 

0 

1.50 

3.00 

4.50 

E 

0 

6.00 

7.00 

8.00 
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Table  4.8 

Measured  Fragment  Velocities 
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Strobe  System  Velocities  l 
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APPENDIX  A 


D.  Lehto 

Calculations  of  Blast  from  Exploding  Pressure  Vessels 


Figure  A-l  summarizes  the  results  of  a  Wundy  Hydro¬ 
code  calculation  for  the  explosion  of  a  6  ft^  spherical  tank  filled 
with  air  at  15°C  and  burst  at  8,000  psig.  The  air  i3  assumed  to  act 
as  an  ideal  gas  with  a  gamma  of  1.4.  Tank  mas 3  is  ignored;  the  con¬ 
tainer  is  assumed  to  disappear  instantaneously  like  an  ideal  shock 
tube  diaphragm.  Table  A-l  summarizes  the  initial  conditions. 

Figure  A-l  compares  the  computed  pressure-distance  re3ult3  for 
the  tank  explosion  with  those  from  a  10.1- lb  TNT  charge,  (if  one 
assumes  that  TNT  energy  is  1018  cals/gm,  then  10.1  lbs  is  the  TNT  energy 
equivalent  based  on  the  Isentropic  expansion  of  the  air  contained  in  a 

6  ft^  air  sphere  compressed  to  8,000  psig  at  15 °C.) 

Figure  A-2  shows  the  space- time  paths  for  the  main  shock,  the 
location  of  the  first  zero  overpressure  (end  of  the  positive  phase), 
and  the  contact  surface  position. 

Figure  A- 3  shows  predicted  pre3sure-tlme  records  at  positions  7-, 

14- ,  and  18-ft  from  the  tank  center.  Contact  surface  arrival  distorts 
the  tail  of  the  7  ft  record. 

Figure  A- 4  presents  the  predicted  pressure-distance  relation  for 

a  1.34  ft^  spherical,  air  filled  tank  burst  at  460  psig.  For  compari¬ 
son,  a  pressure-distance  relation  for  a  0.125  lb  TOT  sphere  is  shown. 
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Table  A-l 

Calculations  for  Exploding  Pressure  Tank 


Volume  (cu  ft) 

Pressure  (psia) 

Gamma 

Wt/Mole 

Tank  Volume  (cu  cm) 

Radius  (cm) 

Compression  Ratio 
Initial  Density  (G/CC) 
Density  in  Tank  (G/CC) 

Mass  in  Tank  (G) 

Internal  Energy  (Erg/G) 

E  Loss  on  Expansion  (Erg/G) 
Total  Energy  Loss  (cal) 

TNT  Equivalence  (Grams) 

TNT  Equivalence  (lb) 


=  6.0 

=  8.0  x  103 
=  1.400 
=  28.964 
=  1.699  x  10s 
=  34.36 
=  5.444  x  102 
=  1.225  x  10* 
=  0.669 
=  1.133  x  105 
=  2.068  x  109 
=  1.726  x  109 
=  4.672  x  106 
=  4.589  x  103 
=  10.12 
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- 10.1  LB  TNT 

- 6  FT3,  8000  RSI 

AIR  SPHERE  (WUNDY  RESULTS) 


O  FROM  SHOCK  TUBE  EQUATION 

.  CURVE  FAIRED  IN  TO  CONNECT 

SHOCK  TUBE  POINT  AND  CALCULATIONS 


SHOCK  FRONT 


_L 

10 


_ L 

20 

TIME  (MSEC) 


FIG.  A-2  SPACE-TIME  PATHS  FOR  EXPLOSION  OF  8000  PSI  AIR  SPHERE 
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A -3  PRESSURE-TIME  RECORDS  FOR  EXPLOSION  OF  8000  PSI  AIR  SPHERE 


OVERPRESSURE  IN  PSI 
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FIG.  A -4  CALCULATED  SHOCK  OVERPRESSURES  FROM  1 .34  FT3  TANK  BURST  AT  460  PSI 
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APPENDIX  B 
by 


R.  A.  Lorenz 


A  Short  Computer  Code  to  Estimate  the  Early-Time  Fragment 
Velocity  and  Position  vs  Time  for  a  High  Pressure  Spherical 

or  Cylindrical  Tank  Rupture 


The  work  summarized  in  this  report  is  the  result  of  the 
need  to  predict  the  fragment  velocity  and  position  vs  time  while 
designing  a  series  of  high  pressure  tank  rupture  experiments. 

The  Spherical  Model 

The  most  important  faotors  contributing  to  the  motion  of  the 
fragments  are: 

(a)  the  driving  pressure  immediately  behind  each  fragment 

(b)  the  jetting  of  compressed  gas  between  the  fragments 
which  results  in  a  lower  driving  pressure,  and 

(c)  the  accumulation  of  compressed  air  ahead  of  the  frag¬ 
ments  and  in  the  jets  which  absorb  energy  and  tend  to 
slow  the  fragments  down. 
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In  constructing  a  model  to  account  for  these  factors  and  to 
keep  the  calculations  as  simple  as  possible,  the  following  basic 
assumptions  were  made: 

(a)  The  total  area  of  all  the  fragments  remains  constant 
while  their  volume  is  negligible. 

(b)  Air  turbulence  and  fragment  tumbling  are  forbidden. 

(c)  At  each  instant,  the  pressure,  internal  energy,  and 
density  of  the  compressed  gas  is  uniform  throughout  the 
volume  of  gas. 

(d)  The  compressed  gas  Inside  the  tank  as-  well  as  the 
surrounding  medium  is  a  const  ant -gamma  ider»l  gas. 

The  geometry  used  in  the  calculation  (Fig.B-$  follows  from 

assumption  fa)  above.  Two  regions  are  defined,  one  for  the  solid  angle 
subtending  the  fragments,  anef  the  other  for  the  jetting  gas.  As  tne 
fragments  move  outward  in  space,  a  larger  fraction  of  the  total 
solid  angle  becomes  available  to  the  jetting  gas  since  the  total 
area  of  the  fragments  remains  constant.  Note  also  that  the  outer 
surface  of  the  jetting  region  is  constrained  to  lie  on  a  sphere. 


The  volume  of  the  system  is  determined  as  follows,  using 
the  notation  of  Figure  B-l: 

Volume  of  fragment  region  =  (-2. ^  nr.^  (l) 

'rs  f 


Volume  of  jetting  region  =  ^1 - (2) 

where  rG  is  the  initial  radius  of  the  tank.  These  calculations  are 
found  in  lines  420-425  of  the  program  on  page  B-10.  (See  also  lines 
390  and  '410-415.) 

As  a  result  of  assumption  (b)  above,  the  motion  of  the  system 
is  restricted  to  radially  outward  motion.  Therefore,  the  key 
velocities  are:  vg,  the  velocity  of  the  fragments;  v^,  the  velocity 
of  the  outer  hour*! ary  of  the  fragment  region;  and  v2,  the  velocity 
of  the  outer  boundary  of  the  jetting  region. 

The  uniform  density  of  assumption  (a)  above,  results  in  a  simple 
expression  for  the  velocity  of  the  gas  at  any  radius  as  a  function 
of  the  velocity  at  the  outer  boundary.  Consider  first  the  time 
rate  of  change  of  the  ''o.  me  of  an  expanding  sphere  in  general: 


dV 

Wc 


4nr2v  =  7  (-|  TTr3)v  =  ~  Vv 


(3) 
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where  v  Is  the  velocity  of  the  outer  boundary  of  the  sphere.  Now 
consider  the  time  rate  of  change  of  the  density  in  that  sphere: 


do  m  dV 
HT  =  "  ^2  37 


where  use  has  been  made  of  Equation  (3).  But  since  o  must  be  constant 
throughout  the  volume,  then 


1  dp  3v 

o  r 


must  also  have  the  same  value  for  any  smaller  sphere  of  constant 
mass  within  the  original  sphere.  Therefore,  it  follows  that 

3v(x)  _  3v(r) 

x  r 


v(x)  =  (£)  v(r)  (4) 

The  kinetic  energy  of  the  gas  in  either  region  of  Figure  (B-l) 
can  now  be  calculated  by  making  use  of  Equation  (4).  First  an 
expression  for  the  spatial  rate  of  change  of  mass  inside  a  constant- 
density  sphere  is  needed: 


dm  _  dV  _ 
dx  ~  D3x  " 


4ttx<: 


(5) 


Then  the  total  kinetic  energy,  T,  of  the  gas  is  found  to  be: 

T  -  f  i  v(x)2dx 

o 

-  fr  i  (o  4nx2}  (p  v2(r)j  dx 

=  fV  Ax  =  Z'.pfW  (pj 

=  r3v2(r)  =  (4  nr3)  ov^ (r) 

T  =  JG  inv2(r)  (6) 
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where  use  has  been  made  of  Equations  (4)  and  (5).  Note  that  the 
result  depends  only  on  the  boundary  velocity  and.  total  mass  within 
the  region.  The  kinetic  energy  is  calculated  in  lines  4*15-455  of  the 
program. 

The  uniform  pressure  and  internal  energy  assumption  (c) 
together  with  the  ideal  gas  assumption  (d)  allows  the  pressure 
to  be  calculated  from  the  expression  for  the  internal  energy  of  an 
ideal  gas,  U,  in  terms  of  Pressure,  P,  and  Volume,  V, 

v  =  a.  (7) 

y-l 

when  the  internal  energy  and  volume  are  known.  Equation  (7)  is  also 
used  to  calculate  the  initial  energy  U0  of  the  system.  The  instan¬ 
taneous  internal  energy  of  the  system  at  any  time  is  defined  simply 
as 

U  m  uo  -  (Total  kinetic  energy)  .  (8) 

Equation  (7)  is  employed  in  lines  195  and  465  of  the  program,  while 
Equation  (8)  is  used  in  Line  460. 

Equations  of  Motion 

In  deriving  the  equations  of  motion,  it  is  necessary  to 
know  the  amount  of  mass  incorporated  into  the  system  due  to  the 
outer  boundaries  of  the  regions  sweeping  up  the  ambient  gas  before 
them.  This  change  in  mass  is  given  by 


Am  =  pfl  AV  =  o*  (4nr  ;  Ar 


or:  Am  »  4tt  q^t  vAt 


(9a) 

(9b) 


where  m  is  the  mass  of  the  gas  in  the  region  and  r  is  the  radius  of 
the  region's  boundary.  Equation  (9a)  is  used  In  lines  335  and  375  In 
the  program. 

Considering  the  jetting  region  first,  let  P  be  the 
pressure  of  the  compressed  gas  and  P,  be  the  pressure  of  the  ambient 
gas.  Then  the  equation  of  motion  for  this  region  becomes: 

(mv)  =  (P-PA)  •  4irr2 


mAv  +  vAm 


=  (P-Pfi)  *  4rrr^ 


or  mAv  =  4TTr2  £(P-PA)  -  pav2J  At 


(10) 
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where  use  has  been  made  of  Equation  (9b).  The  left-hand  side  of  the 
equation  is  the  total  change  in  momentum  of  the  gas  inside  the  region. 
Employing  Equations  (4)  and  (5),  this  term  becomes: 

=  f  af  C(vf(x)  -  v1(x)]  dx  *  f  o  (4ttx2)  [|vf(r)  -  -v^r)]  dx 
o  o 

=  P  [vf(r)  -  vi(r)]  J  x3dx  =  p  Av  (|-) 

o 

=  3  (4  nr3^  o  Av 

mAv  -  mAv  (ll) 

where  Av  is  the  change  in  boundary  velocity.  It  should  be  noted  that 
the  density  o  in  the  above  derivation  contains  an  appropriate  constant 
factor  to  account  for  the  fact  that  only  a  fraction  of  a  total  sphere 
Is  to  be  Integrated  over  so  that  the  ra  in  Equation  (11)  is  truly  the 
total  mass  of  compressed  gas  in  the  region.  Using  Equation  (11), 
Equation  (10)  finally  becomes: 

3  m&v  =  4irr2  [(P-Pa)  -  pa  v2]  4t 

. ^2r(P-PA).,Av2]at 


This  calculation  can  be  found  in  line  350  of  the  program. 

A  similar  derivation  can  be  made  for  the  equation  of 
motion  for  the  fragment  section  except  that  the  equation  corresponding 
to  Equation  (H)  must  also  include  a  term  for  the  change  in  momentum 
of  the  fragments  themselves.  Thus  the  equation  corresponding  to 
Equation  (12)  becomes: 

|  mAv  +  Ms  Av  =  4TTr2  [(P-PA)  -  nAv2  ]  At 


(12) 

(13) 


where  Mg  and  r  are  the  mass  and  radius  of  the  fragments.  This 
calculation  is  made  in  lines  300-305  of  the  program. 

The  change  in  radius  is  given  by  the  formula 

r  =  r  +  vt  +  4vt2 
o 
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or  Ar  =  (v  +  ^Av)  At  (15) 

This  calculation  appears  in  liner  310  and  of  the  program. 
Corrections 

In  order  to  conserve  macs  and  energy,  a  few  small 
corrections  must  he  made  to  the  above  calculations.  First  consider 
the  ambient  gas  which  is  compressed  and  incorporated  into  the  system 
as  the  regions  expand  outward.  The  mass  of  the  added  gas  ha3  already 
been  given  by  Equation  (9a).  If  the  gas  is  compressed  adiabatically 
from  ambient  pressure  Pa  to  the  pressure  P  of  the  compressed  gas,  then 


(16) 


wb^re  f  is  the  fraction  of  solid  angle  occupied  by  the  region,  and 
Ar  is  the  increase  in  radius  due  to  the  newly  compressed  gas.  The 
results  of  Equations  (15)  and  (16)  must  be  summed  to  obtain  the 
position  of  the  region* s  new  boundary.  Equation  (16)  is  used  in 
lines  330  and  370  of  the  program. 

Since  the  velocity  is  computed  by  Equation  (13)  or  (14) 
for  the  old  boundary,  the  velocity  of  the  new  boundary  must  be  found 
using  Equation  (/J).  Therefore,  the  velocitjr  of  the  new  boundary 
will  be  somewhat  greater  than  the  value  resulting  from  Equation  (13) 
or  (1*0.  This  calculation  is  found  in  lines  325  and  365  of  the 
program. 


An  additional  quantity  of  gas  is  being  incorporated  into 
the  system  simply  because  the  Jetting  region  is  expanding  sideways 
as  well  as  outward.  The  mass  concerned  is  represented  by  the 
shaded  area  in  Figure  (B-2).  The  added  volume  can  be  thought  of  as 
the  volume  of  the  shell  between  rO  and  ro’  multiplied  by  the 
fractional  change  Af  in  solid  angle  of  either  the  jetting  or  fragment 
region.  Then  the  additional  mass  is  given  by 
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This  calculation  is  included  in  line  395  of  the  program. 


The  final  correction  considered  is  the  addition  of  the 
Internal  energy  of  the  Incorporated  ambient  gas  to  the  total  energy 
of  the  system.  Using  the  results  of  Equations  (9a)  and  (17)  for 
the  acquired  mass  Am,  the  Increase  in  internal  energy  becomes 


AU 


(‘V)  -  f. x  (|=) 


(18) 


where  use  has  been  made  of  Equation  (7).  This  calculations  is 
performed  in  line  405  of  the  program. 

Cylindrical  Equations 

For  the  cylindrical  calculations,  the  same  model  is  adopted 
as  for  the  spherical  calculations  except  that  axial  symmetry  is  to 
be  used.  Figures  B-l  and  B-2  remain  valid  for  the  new  symmetry .  The 
essential  differences  are  that  the  volume  is  now  nr2h  instead  of 
^nr^  and  the  surface  area  is  2nrh  instead  of  4-rrr2,  where  h  is  the 
height  of  the  right  circular  cylinder. 

The  derivation  of  the  cylindrical  equations  is  similar 
to  that  of  the  spherical  equations,  so  that  only  the  results  will 
be  listed  below  with  the  equation  numbers  100  greater  than  the 
corresponding  spherical  equation  numbers.  The  position  of  the 
1  calculations  in  the  cylindrical  program  is  identical  to  their 
counterparts  in  the  spherical  program.  An  asterisk  (*)  following  an 
equation  number  indicates  that  the  form  of  the  equation  is  the  same 
for  both  symmetries. 


Volume  of  fragment  region  = 

ft) -a 

(101) 

Volume  of  jetting  region  =  | 

(i  -  )  "r2Sh 

(102) 

g  .2  Vy 

cT5  r  vv 

(103) 

< 

'x' 

n 

< 

(104) 

37  =  o2TTXh 

(105) 

p 

T  =  r-rav  (r) 

(106) 
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(107)* 

U  =  TJ  -  (Total  Kinetic  Energy) 

(108)* 

Am  =  oA  (2rrrh)  Ar 

(109a) 

Am  =  (2^rh)  vAt 

(109b) 

mhv  =  2rrrh  [(P-PA)  -  pAv2  ]  At 

(no) 

^  p 

mAv  =  mAv 

(111) 

mAv  =  2nrh  [(P-PA)  -  oAv2]  At 

(112) 

Av  .  3trh  [(P-PA)  -  »AvS]4t 
m 

(113) 

3Ttrh  [(P-PA)  -  pAv2]At 

(114) 

(m  +  1  r~  Ms) 

Ar  =  (v  +  fAv)  At 

(115)* 

,  /MVy 

Ar  =  Ar  Up-I 

(116)* 

Am  =  pATih  (r22  -  r12)  Af 

(117) 

» ■  (“) 

(118)* 

No  account  has  been  taken  of  the  motion  of  the  cylinder 
ends  or  of  the  pressure  release  through  them.  Therefore,  the 
calculations  apply  only  to  the  central  portion  of  the  cylindrical 

tank  where  end-effects  can  be  ignored.  With  this  restriction  in 
mind,  h  can  arbitrarily  be  set  equal  to  unity  in  the  calculations 
since  the  relevant  quantities  become  mass/length  and  volume/length. 
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Accuracy 


The  accuracy  of  the  calculations  will  ultimately  have  to 
be  determined  by  comparison  with  the  experiments  which  have  not  yet 
been  performed.  All  that  can  be  said  a  priori  is  that  when  the 
driving  pressure  is  much  greater  than  ambient'  pressure,  the  defects 
in  the  model  are  expected  to  be  negligible  and  therefore  the  cal¬ 
culated  results  (position  and  velocity  vs  time)  should  be  good.  As 
would  be  expected,  the  driving  pressure  becomes  comparable  to  ambient 
pressure  very  soon  after  the  fragment  velocity  reaches  its  maximum 
value.  By  this  time,  essentially  all  of  the  energy  of  the  system 
has  been  transformed  into  the  kinetic  energy  of  the  fragments  and 
compressed  gas.  Therefore,  as  long  as  the  fragment  velocity  is  high, 
the  momentum  of  the  fragments  will  carry  them  along  and  the  defects 
in  the  model  can  again  be  expected  to  be  negligible.  It  is  when 
the  fragment  velocity  begins  to  drop  significantly  below  the  maximum 
value  that  the  accumulated  effects  of  the  model  defects  should  become 
noticeable. 


The  Programs 


The  programs  are  written  in  FORTRAN  IV  language  for  use  on 
the  CDC  6400  computer  through  the  remote  teletype  units  at  NOL.  The 
listings  for  the  spherical  and  cylindrical  programs  are  well  annotated 
and  are  given  in  Appendices  B-l  and  B-2,  respectively  and  are  given 
starting  on  pages  B-10  and  B-13  respectively. 

Input  data  are  entered  Into  the  program  between  lines  104 — 155. 

Each  program  contains  running  checks  in  the  pressure  and 
internal  energy  of  the  system  and  will  terminate  execution  if  the 
specified  bounds  are  exceeded.  The  tests  are  made  In  lines  490 — 495 
and  lines  760 — 770  of  each  program. 

In  order  to  enable  a  user  to  abort  a  run  and  return  to  the 
SETUP  mode  of  INTERCOM,  the  time-sharing  system  at  NOL,  a  dummy 
READ  statement  is  executed  at  lines  901,  906,  and  911  of  each  program. 
Any  alphanumeric  character  or  set  of  characters  may  be  entered  after 
which  the  program  will  terminate  normally. 

Results 


The  calculated  position  and  velocity  vs  time,  as  well  as 
the  input  data,  for  the  three  pressure  tanks  are  given  in  Figures  B~3 
through  B~5.  The  burst  pressure  is  listed  on  each  figure,  to  which  ambi- 
lent  pressure  must  be  added  before  entering  it  into  the  program  at  line 
115. 

The  same  value  of  y  =  1.4  was  used  for  the  compressed  and 
ambient  gas  in  all  three  calculations.  If  the  gas  considered  were 
air  rather  than  ideal,  the  real  value  of  y  for  even  the  most  compressed 
initial  condition  used  would  not  be  le9s  than  y  ^  **  1.36. 
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Program  Listing  for  Spherical  Calculation 


LilQ 
10‘S 
1  1  0 
1  15 
100 
IP'S 
IP'S 
1  97 
130 
135 
137 
140 
MS 
150 
155 
1  60 
165 
170 
17P 
175 
177 
1«0 
1R5 
190 
1  95 
POO 
PO  5 
91  0 
PI  5 
PPO 
PP5 
P30 
93  5 
940 
PA5 
947 
950 
95  5 
P60 
965 
97  0 
97  5 

oo( 3 


90=13.5  ^ _ 

PO =90*P.54  ^ — 
PO =801 4.7 
P0=P0/1  4.5P.-6 


EJiAG  £J  --TAPF  5—1  NP  itT-i 

4 - Initial  radius 


I! 

!» 

.11 


| inches) 


Initializati 


pressure 

i! 


S‘yASS=14  0.<- 


psi) 
dynes/cm  ) 


Mass  of  Tank  (lbs, 

5MA5s=smass*453.5904  Mass  of  Tank  (grams) 

5  X  A  3  4  3  =  4  .  *  S  M  A  5  5  /  3  . 

i)T=i.-:-6  4 - Timestep  (sec) 

0TP9Mr  =  .05=: -3  < —  Print  interval  (sec) 
tf:.vp  =  3.h:-3  4 - End  of  problem  time  (sec) 

GAMMA  =  1  .4 
GAM  1  =1  .  /  GAM  I"  A 
9H0A  hB=  l  .  PP5»:  -3 
PA  =  1  .01  3P5P6  _ 


PI  =3. 1  4  1  5  9965  35  897  9 
P 1 4  =4 . *PI 
Pl43=Pl4/3 . 

PI  1 63=16. *P I / 3 • 

RH0P4 =9H0AMB*Pl4 
9H3P4  3  =  9H3P4/3 • 

PAG9H0  =PA/CGAt'KA-l  .  )/RHOAMB 

70  =  9  143+90**3  4 - Initial 

A.v  AS  S  0  =R  HO  A  f>  B  *  VO  *CPO/PA)*  *GA  M 1  < 

( 1=  I  JO  =P0 *70/ ( GAMMA -1  .  )  * - 

RS=91  =99=90  <- - 

!?50P=9S*9S 

7FLs=v-:li  =vflp=o. 


Ambient  Density  (gm/cc)  s 
Ambient  Pressure  (dynes/cm 


volume  (00) 

— Initial  Mass  of 

-  n  internal 

-  "  (area  of 


gas  in  tank  (gm) 
energy  (ergs^ 


tank^/4n  (cm2) 


F 1  =  1  . 

F  9  =0  •  4~ 


AMASS=AMASSO 
A  1  MASS  =AMA  SS 
A9MASS  =0 • 

T  =0  .  ^ - 

TPRMT=  -Df/P  •  <- 


Fractlon  of  sphere  in  fragment  region 
"  "  "  "  jetting  " 


Time  (sec) 

-  Printout  time  (sec) 


tfmo  =t-:  \'d+tprm  r 
p=po 
v  =vo 

PRINT  8 8 >  PO  > PA  >  RO  * 70 > GAMMA , 9H0A MB , UO  >PAGRH0 , SM AS S *  A M ASSO *0 T 
8  8  F09MAK/4H  P0  =  1  PO  1  0 . 3^4X3HPA=F  1  0 . 3  >4 X3HR0  =F  1  0 .3 *  4X3H70  =9  10-3/ 

+  7H  GAMMA =0PF5  .P  ,4X5HRH0A  =1  P~  1  () .  3  jA  X3H  LO  =E  1  0  •  3, 4  X7HPAG9H0  =E  1  0  •  3/ 
+  7  H  SMA  SS  =€  1  0  •  3  ^  4  X7  HAM  A  S.SO  =F  1  0  •  3  >4  X3H0T  =£ 1 0. 3/ > 

GO  TO  POO 
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100  DVELl  =Fl*PI163*(Rl**2)*(<P-PA)-RH0Af,'b*\/£Ll**2>*OT 
DVELl  =DVEL1  /CA1MASS+CRS/R1  )*SFAS43> 

dpi  =(v- li+oveli/2. )*ot  Eq.  of  motion  calcs 

ps =rs+dri *crs/ri )  Fragment  Region 

PAPGl  =  CPA/P)**GAM  ‘ 

VEL 1 =(VEL1 +DVEL1 )*C  t .  +  ( PAPGl *DR  1  >/ (R 1 +DR1  )) 

DPI =0R ! * C 1 .+PAPG1 ) 

DAI  PAS  =  F1  *RH0P4*CR1  **2>*DR1 
PI =R 1 +DR1 

VELS=VEL1  *<R$/R1  )  .  . . . 

0  VE  L2  =  (  F2*PI  1  63*  (  R2**2>  *(  <  P  -PA  )  -RH0ANB* VE  L2**2 ) *DT  )/A2NASS 
IFCA2MASS.EQ.0  .  )  DVEL2.=DVELl 

DR2=(VEL2+DVEL2/2. >*DT  w  r  ..  . 

VE  L2  =  C  VEL2+0VE  L2 ) * (  1  •  +  ( PAPG 1  *0 R2  >/ ( R2+D R2 )  )  “O’1'1®11  caJ.CS 

oro=!)rc5+c  i  .  +  papgi  )  Jetting  Region _ 

DA2KAS=F2*RH0P4*(R2+*2  )*DR2 
...S2=R2  +  D  R2 . „ 

Fl N=RS02/RS**2  . 

on  =rhop43* crp**3 -ri  **3 ) *cfi -fi n >+o a  1  |vAs+DA2tiAs  l  Increment  total 


AFASS  =AFASS+DK 

IJO  =  U0_+  PAGR HO  *0  N _  _  _  _  ___ 

FI  =F1N 
F2  =  l  •  -Fl 

V  0  LI  =PI43*(R1 **3 )  *F  1 
VO L2 =P 143* ( R2**3  > *F2 
VN=V0L1  +V0L2 

A  1  t'ASS=AFASS*<V0Ll/VN) 

A2  FA  S  S  =  A  frA  S  S  -A  1  FA  S  S  _  _ 

n =.3*A1 FASS+VEL1  **2 
T2=.3*A2MASS*VEL2**2 
TS  =  .  5  ♦SN'AS  S*VE  LS**2 
UM  =  UO  —  C  T  1  <-T2.+TS) 

P=P*(V/VN)*(1JN/'J ) 

V  =  VN 

U^UN_  _  _ _ _  _ 

T=T  +  DT 

I  F  (  r  «  GE •TPRNT )  G0  T0  200 
IFCP.LE.PA.OR.P.Gr.PO)  G0  TO  200 
IFCU.LE.O.)  GO  TO  200 

CTTT0  100  _  _  “ 

200  A  A  =  A  N  A  S  S  /  A  VA  S  SO" 

PRINT  ROW  T  >  RS  *  VELS  jP>  V»  AA 
PRINT  802  >  F1,R1,VEL1  / 

PRINT  803>  F2>R2,VEL2  C 

IFCT. GE .TEND)  GO  TO  700  J 

IF(P.LE.PA.OR.P.GT.PO)  GO  TO  710 

IFCU. LE.O.)  GO  TO  720 
TPRNT=TPRNT  +  DTPR.NT 

GO  TO  100 


ma3s  and  energy 


Reapportion  the  mass  between 
the  two  regions 


Calculate  new  pressure  and 
Internal  energy 


Test  for  Printout 


Parameters  are 
described  on  page  B-3 


Printout  routine 
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900  700  PRINT  808 

901  READ  (5*807)  AA 

902  STOP  77 

90S  710  PRINT  808 

90  6  READ  (5*807)  AA  < - — ^Dnmmy  rea 

907  STOP  1  1 

910  720  PRINT  808 

911  READ  (5*807)  AA  WT 

912  STOP  22 

92  0  801  FOR NAT (/I X 1  P  5  E 1 3*5  *  0  P  F  6  •  3) 

93  0  80  2  FORMAT (4  X3HF1 =F6 • 4  *  1 X 1 P2F 1 3 • 5 ) 
940  803  FORMAT (4X3HF2=F6 *4  *  1 XI P2E 1 3 . 5 ) 
950  807  FORMAT ( A 1 ) 

960  808  FORMAT (//3H  A=) 

999  END 


•Dummy  read  for  abort  purposes 


PRINTED  OUTPUT 


LINE  720 


VELS 


Time  (sec) 

Radius  of  Fragments  (cm) 

Velocity  of  Fragments  (cm/sec) 

Pressure  of  Compressed  Gas  (dynes/cm2) 
Volume  ”  ”  "  (cc) 

Ratio  of  Mass  of  Gas  to  Original  Mass 


LINE  730 


FI  =  Fraction  of  Solid  Angle  for  Fragments 
Rl  *  Outer  Radius  of  Gas  in  Fragment  Region  (cm) 
VEL1  =  Velocity  at  Boundary  of  Fragment  Region  (cm/se< 


LINE  7^0 


F2  « 
R2  = 
VEL2  = 


Same  as  LINE  730  but  for 
the  Jetting  Region 
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Program  Listing  for  Cylindrical  Calculation 


JOG 
1 05 
110 
115 
120 
125 
125 
127 
130 
135 
137 
140 
145 
150 
155 
150 
155 
170 
175 
177 
1  80 
185 
190 
1  95 
200 
205 
21  0 
21  5 
220 
295 
23  0 

23  5 

24  0 
945 
94  7 
950 
255 
960 
965 
970 
975 
9P’0 


_ FRAG(  I.NPur* OUTPUT* TAPE 5 

26=5.25’  ^—-Initial  radius 


=  input ) 


If 


Initialization 


| inches) 
psi) 

6< — it  ti  (dynes/cm  } 

Mass  of  Tank/Length  (lbs/in) 

< - Mass  of  Tank/Length  (grams/cm) 


-Initial 
« 


pressure 

I! 


20  =R0*2 • 54 < 

PO  =4  74 .7^ _ 

P0=P0/1  4 . 5 F 
SMASS=.31{— 

SFASS=S  h ASS *45 3. 5 924/2. 54 
SFAS32=3. *SMASS/2- 

dt=i.e-6^ - Timestep  (sec) 

o  t  p  rn t  = .  o  5 E  - 3  - Print  Interval 


-End-of -Problem 


TEND  =3  -E  -3^ _ 

GAFFA  =1  .4 
GAM1  =1  ./GAFMA 

RH0AMB  =1  . 22 5E -3  f - 

PA  =  1  .01  325E  6  - - 

PI =3. 14159265358979 

P 1 2  =p  l  *2 . 

PI3=PI*3. 

RH0P=RH0AFB*PI 

RH0P2=RH0P*2« 

PAGRHO  =PA/  (  GAFiFA  -1  .  J/RH0AFB 
00=PI*R0**2  { - 


(sec) 

Time  (sec) 


-Ambient  Density  (gm/cc) 

—  "  Pressure  (dynes/cm  ) 


Initial 

_ n 


AFASSO=RH0AF8*VO*(PO/PA  >**GAF1 

U=1JO=PO*VO/C  GAFMA-  1  •  )  ^ - 

RS=R1 =R2=R0 
RSO =RS 

V-'  LS  =VE  L 1  =VFL2  =0. 

f  i  =  i .  (. - Fraction  of  Cylinder  in  Fragment 

F  2  =0.  6 - 


Volume  (cc) 

Mass  of  Gas  in  Tank  (gm] 
Internal  Energy  (ergs) 


it 


ti 


ii 


it 


Jetting 


region 

it 


A  PASS  =A  FASSO 
A  1  t1  ASS  =AMA$S 
A2FASS  =0. 

T  =0 .  < - 


-Time  (sec) 
- Printout 


time  (sec) 


TPRNT  =  -DT /2.  <; - 

TE NO =TEND+ TPRNT 
P=P0 
V  =90 

PRINT  8  8*  PO*PA*RO*VO*GAFMA,RH0AFB*  (JO  *  PAGRHO  *  SPASS  >  AMASS0  >0T 
88  FORMAT C/4H  PO  =1  PE  1  0 . 3  *4 X3HPA  =E  1  0 . 3  *4  X3HR0  =E  1  0 . 3  *  4 X3HV0  =E  1  0 . 3/ 

+  7H  GAFFA  =QPF5 . 2* 4X5HRH0A  =  1  PE  1  0  •  3  *4  X3HU0  =E  1  0  • 3  *4  X7HPAGRHQ=E  10*3/ 

+  7 H  SMASS=E10.3*4X7HAFASS0=El0.3>4X3HDr=£10.3/> 

GO  TO  200 
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^9  *00  D^Lt  =Fl  *PI3*Rl  *(<P-PA> -RH0AtoB*VELl  **2>*DT 

3.0  dr! --wTu iiwJi J5?n£S/R ' >,SMAS32 ’  fragment°reelonal°B 

315  RS  =  RS+0R1  +CRS/R1  )  Sl*.8Vent. , .£?$!?** 

320  PAPG1  =<PA/P)**GAN1 

^  nn^no^/ 1+DVEL1  >*<1  *+{pAPGl*DRl  )/(Rl+0Rl  )) 

330  DR1 =DR1 *( 1 ,+pAPGl ) 

335  DAI FAS=F 1 *RH0P2*R1 *DRl 
340  R1  =R1  +DR1 
345  VE LS =VE L 1 *< RS/R 1  ) 

350  DVEL2  =  <F2*Pi3*R2*UP-PA>-RH0AI*B*VELp**2}*DT>/A2MAe*<5 
355  IFCA2MASS.EQ.OO  DVEL2=DVEL  1  ur>/A2MA.S 

360  DR2= ( VE L2+D VEL2/2 • )*DT  Eq,  Of  motion  CalCS 

365  VE L2  =  ( VE L2+DVE L2  )  *  (  1  *+(PAPGl  *D R2  ) / ( R2+ 

370  DR2=0R2*C  1  .  +  PAPG1  )  . 

375  DA2NAS =F2*RH0P2*R2*DR2 
380  R2=R2+DR2 
390  F1N=RS0/RS 

400  a^ass=anass+dm ~^***^*^*  ~fi n  >+d a  1  MAs+o A2KAS  / Increment  total  mass 

405  liO.=  UO+PAGRH0+Dy, _  \  30(3  enerS7 

410  F1=F1N  “~Y —  —  — * 

S  I  rxizz mass  between 

430  VN=V0Ll  +  V0L2  (  regions 

435  A 1  MASS  =A MASS* (VOL  1 /VN )  \ 

440  A2^S  =ANASS_-A1  t^ASS  _  J 

445  T 1 =*25*A1 MASS*VELl **2  1  \ 

450  T2 =.25*A2NASS*V£L2**2  j 

tl  S^r^rsf*2  i  «“e  —  «™*  - 

465  P=P*(V/VN>*<  (JN/U)  (  ** 

470  V=VN  \ 

475  U=UN _  ) 

480  T  =T +DT  ^ 

4on  rp-fi*Pr‘lARI^T)  60  T0  200  sTest  for  printout 

490  I FCP . LE • PA. 0R. P. GT . PO )  G0  T0  200  \ 

495  IFCU.LE.Q. )  G0  TO  200  ) 

500  G0  P0  100  —  - / 

700  200PAA=AMASS/AMASS0 

720  print  801  *  i  > rs>vels*p>\/,aa  Printout  routine 

730  print  802 j  fi,ri,veli  >  Parameters  are  “ 

740  print  803>  F2,r2,vel2  \  described  on  page  B-3 

750  IF(T.GE.TEND)  G0  T0  700  ^ 

760  IF( P . LE • PA • OR. P . GT  •  PO  )  G0  T0  710 
770  IF<U»LE.)  G0  T0  720 
780  TPRNT  =TPRNT+DTPRNT 
790  G0  T 0  5  00 
900  700  PRINT  808 


Test  for  printout 


Printout  routine 

Parameters  are 
described  on  page  B-3 
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901 

902 

905 

906 

907 
910 
91  1 
912 
920 
930 
940 
950 
960 
999 


READ  (5*807)  AA 
STOP  77 
710  PRINT  808 
REhD  (5*807)  AA 
STOP  1 1 
720  PRINT  808 
READ  (5*807)  AA 
STOP  22 

801  F0RNAT(/1X1P5E13.5*0PF6.3) 


Program  Windup 
Dummy  read  for  abort  purposes 


802 

803 

807 

808 
END 


FORMAT  (4X3HF1  =F6« 
FORMAT (4X3HF2  =F6< 
FORMAT ( A1 ) 

FORMAT (//3H  A=) 


4*1X1 P2E 13*5 ) 
4*1X1 P2E 13*5) 
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■  =  RADIUS  OF  TANK  FRAGMENTS 

■,  =  RADIUS  OF  COMPRESSED  GAS  AHEAD  OF  FRAGMENTS 

■2  =  RADIUS  OF  JETTING  GAS 


FIG.  B-l.  BASIC  GEOMETRY  FOR  THE  CALCULATIONS 
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O  CN  O  CO  O 

CO  00  00  CO  CM  CN 


CN  O  CO  sO  ^  CN  O 
CN  CN  - —  . —  »—  *— 


(03S/UJ0)  e0i  x  A 


FIG.  B-3.  FRAGMENT  RADIUS  AND  VELOCITY  VS  TIME 
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-<r  cn 

oo  oo 


©  00  O  Tf 

CO  CM  CM  CN 


CN  °  2  2  2  2  2 


(oas/tuo)  ^0|X  /\ 
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B-5.  FRAGMENT  RADIUS  AND  VELOCITY  VS  TIME 
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